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ABSTRACT 

A  new  method  for  absolute  calibration  of  electroacoustic  trans- 
ducers is  described  which  permits  low  frequency  reciprocity  calibration 
in  a  test  apparatus  of  small  dimensions.  The  method  takes  advantage  of 
the  reduced  wave  speed  in  a  cylindrical  water  column  bounded  by  a 
compliant  PVC  tube  to  reduce  the  length  of  the  standing  wave  to 
dimensions   which    can    be   easily   handled    in   the   laboratory. 

-  Results  of  a  limited  series  of  experiments  show  a  mean  reproduci- 
bility of  ±0.8  decibels  in  a  one-meter  long  tube  at  frequencies  between 
750  and  1100  Hz  and  sound  speeds  from  330  to  365  meters/second.  The 
hydrophone  voltage  sensitivities  obtained  by  the  method,  however, 
systematically  differ  by  1.5  decibels  from  those  measured  by  comparison 
with  a  standard  hydrophone;  additional  research  will  be  required  to 
resolve  this   discrepancy. 
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I.    WHY   COMPLIANT  TUBE    RECIPROCITY  WORKS 

A.      ACOUSTIC  CALIBRATION   AND   RECIPROCITY 

There  are  two  general  classes  of  underwater  acoustic  transducer 
calibration:  absolute,  or  primary,  calibration,  and  comparison,  or 
secondary,  calibration.  The  latter  presumes  the  availability  of  a  sound 
field  of  known  amplitude  measured  by  a  calibrated  transducer,  but 
absolute  calibration  can  be  performed  using  only  voltage  measurements 
and  primary  measurements  of  the  test  system  environment  in  a  sound 
field  of  unknown  amplitude.  Once  the  voltage  sensitivity  of  a  standard 
hydrophone  has  been  established  in  a  particular  frequency  range  and 
environment  by  absolute  calibration,  it  may  be  used  as  a  reference  for 
comparison  calibration  of  other  hydrophones  and  of  projectors;  this 
absolute  calibration    is   most  often   done   by   the    reciprocity   method. 

1  .       Fundamentals   of    Reciprocity   Calibration 

The  reciprocity  method  for  absolute  calibration  is  based  on  the 
properties  of  linear,  passive  systems  whose  action  can  be  described  by 
the  two-port  equation  (1.1),  which  relates  the  acoustic  pressure  and 
velocity  (p,u)  at  a  transducer's  mechanical  terminal  and  the  voltage  and 
current  (E,l)  at  its  electrical  terminals  by  a  symmetric  or  antisymmetric 
matrix   of   constant   coefficients;    such    a    system    is   called    reciprocal. 


.±$      ZE 


(1.1) 


Hunt  [Ref.  1]  derives  this  equation  from  first  principles  and  shows  how 
the  matrix  elements  may  be  evaluated  for  common  transducers.  Foldy 
and  Primakoff  show  [Ref.  2]  that  the  short  circuit  transmitting  current 
response  S  and  the  open  circuit  hydrophone  voltage  response  M  of  a 
reciprocal   transducer  T   are    related   by   a   constant  factor 


where  J  is  called  the  reciprocity  parameter  and  is  characteristic  of  the 
test  medium  and  its  boundaries.  Bobber  has  demonstrated  [Ref.  3] 
that  J  is  equivalent  to  the  acoustic  transfer  admittance  between  an 
arbitrarily  defined  transmitting  surface  T  generating  a  volume  velocity 
U  and  an  arbitrarily  defined  receiving  surface  H  within  the  test  envi- 
ronment exposed  to  an  average  pressure  P,  and  exists  as  a  function  of 
the  test  medium  and  its  boundaries  for  any  arbitrary  linear,  passive, 
and    reversible  medium. 

When  two  test  transducers1  are  located  at  points  in  the 
medium  where  the  pressure  amplitudes  are  equal,  the  reciprocity  param- 
eter allows  expression  of  the  signal  voltage  of  a  test  hydrophone  H 
placed  at  one  point  as  a  function  of  the  transmitter  current  response  of 
a  reversible  transducer  T  placed  at  the  other  point  and  driven  to 
create   the   sound   field; 


Eth  =  MHPH  =  MHSTIT  (1  .3) 

where  I  is  the  drive  current  input  to  T.  If  the  sound  field  generated 
by  T  is  then  replaced  by  a  field  which  is  identical  at  the  locations  of  T 
and  H,  but  is  generated  by  an  separate  projector  P,  the  outputs  of  T 
and  H  in  this  new  field  will  be  related  by  equation  (1.4).  Combining 
this  equation  with  (1.3)  produces  the  calibration  equation  (1.5),  which 
allows  the  computation  of  the  test  hydrophone's  voltage  sensitivity  H 
directly  from  voltage  measurements  made  in  the  three  source-to-receiver 
configurations:      T-to-H,    P-to-H,    and    P-to-T. 


:To  avoid  confusion,  the  term  test  transducers  will  always  be  used 
when  referring  to  both  transducers  used  in  the  reciprocity,  calibration. 
The  term  reversible  transducer,  or  the  ANSI  designator  "T",  will  be 
used  to  refer  to  the  reciprocal  transducer  used  both  as  a  hydrophone 
and  as  a  projector  during  the  calibration;  test  hydrophone,  or  "H"  will 
designate  the  transducer  used  only  as  a  hydrophone;  and  projector,  or 
"P"  will  identify  a  third  sound  source  used  to  create  the  sound  field 
used    in   the  test. 
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Li  =  *«£  =  JJx  (1.4) 

EPH      MH°      .  MH 

2         E     E  CI. 5) 

Mh  =  J  ^^ 

EpT  't 


2.      Reciprocity   Test  Geometries 

Most  reciprocity  calibrations  are  made  in  one  of  the  two 
earliest  system  geometries  for  which  the  reciprocity  parameter  was 
computed:  free  field  [Ref.  4],  or  close-coupled  [Ref.  5].  The  tech- 
niques for  carrying  out  these  tests,  and  the  limitations  required  to 
ensure  the  validity  of  the  test  are  described  by  the  ANSI  standard  for 
underwater  calibration  [Ref.  6].  The  free  field  test  geometry  requires 
a  volume  of  water  large  enough  to  permit  spherical  spreading  of  the 
acoustic  waves  undisturbed  by  boundary  reflections.  This  implies  tanks 
of  dimensions  greater  than  ten  wavelengths  (about  15  meters  at  1  kilo- 
hertz)  for  underwater  testing.  Coupler  reciprocity  is  limited  by  the 
requirement  to  encase  the  entire  active  surface  of  two  test  transducers 
in  a  rigid-walled  cell  of  dimensions  less  than  1/10  of  a  wavelength; 
useful  underwater  low  frequency  transducers  limited  to  such  dimensions 
can    be  difficult  to   build. 

The  practice  of  reciprocity  cal.ibration  is  not,  however, 
constrained  to  the  two  conventional  geometries.  Rudnick  [Ref.  7] 
discusses  three  compact  geometries  with  known  reciprocity  factors; 
Bobber's  reference  text  on  measurements  [Ref.  8]  describes  several 
more.  These  unusual  test  environments  extend  the  practice  of  reci- 
procity calibration  beyond  the  limits  of  of  the  two  conventional  geome- 
tries. A  summary  of  these  methods  is  included  in  Table  I;  each  of  the 
reciprocity  parameters  listed  is  the  quotient  of  an  area  characteristic  of 
the  system  geometry,  and  of  the  specific  acoustic  impedance  of  the 
medium. 

The  compliant  tube  test  system  described  in  this  paper  adds 
another  option  to  the  array  already  available  to  the  experimenter.  It 
has       particular      value       to      the       underwater      acoustician       because       it 
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TABLE   I 
Some   Reciprocity   Calibration   Systems 

Geometry  Reciprocity  Parameter  Reference 


Free  field 


2dX 
P0c 


Close— coupled  r 


Plane— wave 
Cylindrical-wave 
Helmholtz  resonator 


P0c 

2A 
P0c 


2LVdX  8 

P0c 

2rrALf  7 

P0c2Q 


Waveguide,  ?n 

resonant  plane  wave  Poc  ^ 


significantly  lowers  the  frequency  limit  for  underwater  reciprocity  cali- 
bration below  the  practical  limits  of  the  conventional  geometries. 
Unlike  rigid-tube  methods,  it  can  be  accomplished  within  a  convenient, 
inexpensive  test  cell,  and  although  still  limited  by  the  requirement  for 
accurate  measurement  of  the  system  quality  factor,  does  escape  resonant 
reciprocity  calibration's  usual  restriction  to  discrete  frequencies.  The 
"virtual  rigid  walls"  created  by  the  standing  wave  pressure  antinodes, 
and  in  which  the  transducers  are  placed  for  the  test,  can  be  adjusted 
to  any  resonant  length  and  frequency  by  simply  raising  or  lowering  the 
water   level. 

Our    development    of    the    compliant    tube    reciprocity    parameter 
must    begin   from    an    understanding   of    wave  motion    in    the  water    column, 
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and    of    the    oscillatory    expansion    of  the    tube's     circumference.        From 

these    equations,    we    will    develop    the  energy    balance    of   the    oscillating 

system;      this      in     turn     will     provide  the      result     for     the      reciprocity 
parameter. 

B.      COMPLIANT  TUBE  WAVE   EQUATION 

1  .      Waves  ]_n   the  Water 

Lighthill's  development  [Ref.  9]  of  the  wave  equation  for  a 
lossless  fluid  moving  coaxially  in  a  lossless  tube  of  varying  cross 
section  includes  the  case  of  the  compliant  elastic  tube  in  which  our 
experiment  is  performed.  The  momentum  and  mass  conservation  equa- 
tions    for  an    inviscid,    ideal   fluid   are 

( du        du\         do 

+  u—    =  — - 


p^+u—    =  -^  (1.6) 

Vat        dz)         dz 

d(pa)        a  (1.7) 

— —  +    j(PQU)      =  ° 
a  t  dz 

where  p  is  the  acoustic  excess  pressure,  u  the  acoustic  coaxial  particle 
velocity,  and  a  the  local  time-varying  cross-sectional  area  of  the  tube, 
which  varies  slightly  about  its  equilibrium  value  A.  When  fluid  entropy 
is  constant,  both  the  local  density  p  and  the  local  cross  section  a  may 
be  expressed  as  functions  of  the  acoustic  pressure.  (The  latter  relation 
is  just  the  inverse  of  the  function  which  relates  pressure  to  cross 
section  variations.)  Discarding  terms  of  second  order  and  higher 
reduces   these   equations   to 


n  ^-  _^E                                                                                                 (1.8) 

H  at  dz 

d(pa)  du                                                                                       H.9) 

at  H    dz 
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Taking  the  spatial  derivative  of  (1.8)  and  the  time  derivative  of  (1.9), 
combining  the  results  with  the  linearized  equation  of  state  p  =  (p  -p0)c2 
and   discarding   terms   above  first  order  gives   a   wave   equation 


afp_   _2a2P 
dz~  a?" 


(1.10) 


where  the  wave   speed   c   is   defined   by 


•~2=      of—  d(pa))    =    p0  (   -  — 

°\p0A  ap    /o       °  \  p0ap 


L  ^5 
a  ap. 


=  p0     K    +     D 


(1.11) 


(The  subscript  "0"  above  implies  evaluation  of  the  derivative  at  the 
undeformed  wall  position  and  under  adiabatic  conditions.)  The  first 
term  of  equation  (1.11)  is  the  adiabatic  compressibility  or  reciprocal 
bulk  modulus  K;  if  there  were  no  dependence  on  cross  section,  this 
would  reduce  to  the  equation  for  the  sound  speed  in  the  bulk  medium. 
If  the  second  term,  called  the  disten-sibility  D,  is  much  larger  than  the 
adiabatic  compressibility,  the  speed  of  sound  waves  in  the  tube  will  be 
determined  by  the  characteristics  of  the  wall  rather  than  by  those  of 
the  bulk  fluid,  and  will  be  much  slower.  This  in  turn  will  reduce  the 
size  of  the    resonant  tube   required   to  operate  at   a   given   frequency. 

The  harmonic  solution  of  equation  (1.10)  with  a  pressure 
release  boundary  at  z=0  and  a  rigid  boundary  at  z=L,  and  for  frequen- 
cies  below   the   cutoffs   for   azimuthal    and    radial    modes,    is 


p(z,t)  =  P  sin(kz)sin(cjt) 


(1.12) 


with    a   corresponding   velocity   field 
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u(z,t)  = cos(kz)  cos(o;t) 


(1.13) 


The     wave    number     k     is     limited     by    the     boundary     conditions     to     the 
discrete  values 


k    = 


2rr  _    7r(n  +  1/2) 


A 


n=0,  1,  2,  3,.... 


(1.14) 


This  sound  field  has  associated  with  it  both  kinetic  and  potential  energy 
density.  The  kinetic  energy  of  axial  water  motion  in  a  thin  radial  cross 
section   of   length   dz   is 


dEu  =  £fAdz 


(1.15) 


The    potential    energy    in    the    same    section     is    found    by    considering    the 
work   done  by   the  fluid   as   the  tube   boundaries   expand: 


dEp  =  - 


<  r<      < 

p  dV  =  -         p  ]  Adz 


0 


1  da 
A  dp. 


a  a     1    dP 

-  A  dz 

PodpJ 


dp 


(1.16) 
(1.17) 


=  (K  -  D)*-    Adz 

where   C    is   a   coordinate  of  the   periodic  motion. 

2.      Motion   of  the  Wall 

The  distensibility  can  be  calculated  from  the  material  proper- 
ties of  the  tube  wall  and  its  geometry.  In  the  case  of  a  cylindrical  tube 
which  responds  elastically  to  small  excess  stresses,  the  circumferential 
tensioning  force  component  due  to  the  pressure  on  a  thin  radial  cross 
section   of   length    dz   is 


15 


.n/2 

p  cos  8  d9  =  4Rdz  p 

-n/2 


F  =  2Rdz 

J    -n/2 

<5o^    F         =2Rp 
A      2hYdz   "  hY 


(1.18) 
(1.19) 


This  tensioning  force  causes  a  strain  s  defined  by  equation  (1.19). 
Where  p  is  small  and  a  varies  only  slightly  from  its  undisturbed  value 
A,    the  distensibility   is   then 


D  = 


1  da 
A  dp_ 


2R 
hY 


(1.20) 


This  stretching  of  the  wall  along  its  circumference  has  the  characteris- 
tics of  a  damped,  stiffness-controlled  driven  oscillator  (neglecting  long- 
itudinal effects).  The  ratio  of  the  driving  pressure  at  any  cross 
section  and  the  circumferential  strain  velocity  is  the  mechanical  impe- 
dance of  the  cross-sectional  ring  treated  as  a  damped  oscillator  of 
negligible  mass  with  elastic  constant  Y  and  damping  coefficient  R(m) 
[Ref.    10:pp.  10-12] 


£  = 


2P  sin(kz)sin(cjt-0) 


cjW 


w  =  n  r;+(y/«)' 


9  =  arctan    


(1.21) 
(1.22) 


The  polyvinyl  chloride  (PVC)  tube  used  in  this  experiment  is  actually 
not  elastic  over  any  finite  range  of  distortion,  so  the  parameter  Y  is 
not  Young's  modulus;  it  is  only  a  linear  approximation  to  the  plastic 
deformation  that  the  PVC  undergoes  in  harmonic  motion.  Within  this 
linearization,  the  potential  energy  stored  in  a  thin  radial  section  of  the 
wall  as  the  tube  is  expanded  by  a  small  amount  from  its  initial  cross 
section    can   be  found   from   the   force-displacement   integral 
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F—ds  =  Dp2 
o     2 


Adz 


(1.23) 


with   the  force  given   by   equation    (1.18). 

C.      COMPLIANT  TUBE   RECIPROCITY   PARAMETER 

The  arguments  used  by  Garrett,  et  aL  ,  [Ref.  11]  to  derive  the 
reciprocity  parameter  J  for  a  rigid  tube  will  also  yield  an  expression  for 
the  reciprocity  parameter  of  the  free  surface,  compliant  wall  tube.  As 
with  the  rigid  tube,  the  solution  depends  on  relating  the  energy 
content  of  the  system  E  and  the  energy  loss  per  cycle  AE  to  the 
acoustic  transfer  admittance  (reciprocity  parameter)  described  by  equa- 
tion (1.2).  When  the  ratio  of  energy  content  to  energy  loss  rate  is 
known  in  terms  of  J  and  of  primary  measurements,  the  resonance 
quality  factor  equation  (1.24)  (which  can  also  be  found  from  primary 
measurements)    can    be   used   to   evaluate  J. 


Q  = 


fh"f, 


AE/2tt 


(1.24) 


1 .      System    Energy   Content   and    Losses 

The  compliant  tube  system  stores  energy  both  in  the  water 
column  and  in  the  tube  wall.  Kinetic  energy  of  motion  of  the  water 
pressurizes  the  water  locally,  which  in  turn  stresses  and  expands  the 
tube  diameter.  The  energy  balance  exchanges  between  the  three 
degrees  of  freedom  are  governed  by  the  work  equations  for  pressuriza- 
tion  (equation  (1.17))  and  for  expansion  (equation  (1.23)),  and  by  the 
kinetic  energy  equation  (1.15).  Summing  these  equations  gives  the 
total    energy   in    a   thin   cross    section 
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dE  =  dEu  +  dEP  +  dE„     =  A  dz  \  ^+  (K  -  D)|-+  Dp 


=  Adz|^+|"2(K  +  D)|=Adz(     2 


2p0c 


(1.25) 
(1.26) 


where  equation  (1.11)  has  been  used  to  replace  K  and  D  with  p0  c2 . 
The  sound  field  in  the  absence  of  attenuation  is  described  by  equations 
(1.12)  and  (1.13);  substitution  gives  a  form  which  can  be  integrated 
over  the   length   of  the  water  column 


rV 


E  = 


•6 


2    \p0c 


cos2(kz)cos2M)  +—  2  sin:(kz)  sin2(cjf)f  dz     (1  .27) 
2p0c  ) 


ALP  - 
4p0c 


(1.28) 


This    is    identical   to  the    result  for  a    rigid    resonant  tube    [Ref.    11]. 

In  the  steady  state,  the  cyclic  energy  loss  equals  the  work 
done  by  the  driving  transducer.  If  this  transducer  is  located  at  a 
pressure  antinode,  pressure  and  velocity  must  be  in  phase  at  the 
transducer  face  (  since  the  impedance  to  flow  at  a  velocity  node  is  infi- 
nite). The  energy,  transferred  to  the  system  is  the  work  done,  which 
can    be  written   as 


3E 

at 


=  0puonds     =PU  sin2(cjf) 


(1.29) 


where    U     is     the    volume    velocity    amplitude    of     the    transducer    surface. 
Integrating   over  one  cycle  gives 


PU 

AE  =  — 

2f 


(1.30) 
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2.      Reciprocity   Parameter 

Our  calculations  above  have  yielded  both  the  system  energy 
content  and  the  system  energy  loss  (actually,  its  throughput  at  steady 
state).      The  energy   equation   for  Q    (1.24)    immediately   produces 


Q  E        PALf 

-  =-   =1TT— 2  (1.31) 

2tt      AE      U  2p0c 


The  ratio  of  volume  velocity  amplitude  to  pressure  amplitude  that 
appears  in  this  result  is  the  acoustic  transfer  admittance  between  pres- 
sure maxima  (equation  (1.2)),  and  must  also  be  the  reciprocity  param- 
eter for  hydrophones  positioned  at  separate  maxima  when  the  column  is 
driven  at  resonance.  Equation  (1.32)  follows;  again,  it  is  the  same  as 
the    reciprocity   parameter  for  the   rigid   tube. 


PoC 


^Q  (1.32) 


The  solution  of  equation  (1.32)  requires  the  experimental  measurement 
of  three  resonance  parameters:  the  resonance  frequency  f;  the  reso- 
nance quality  factor  Q;  and  the  wavelength'  of  the  resonance  standing 
wave,  in  order  to  compute  the  sound  speed  C.  The  accuracy  of  the 
value  computed  for  J  depends  on  our  ability  to  measure  these  parame- 
ters accurately.  The  methods  for  doing  this  will  be  discussed  in  the 
next  chapter. 

In  using  equation  (1.5)  with  the  result  derived  above,  we 
must  ensure  that  the  pressures  and  voltages  implied  are  indeed  the 
maximum  amplitude  values  at  the  antinodes.  The  test  transducers 
output  will  be  proportional  to  the  average  sound  field  over  their  active 
surfaces,  which  at  its  maximum  will  be  less  than  that  which  would  be 
generated  by  a  uniform  field  at  the  maximum  pressure  amplitude  P  by 
the  factor 
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-   1  rf\  i  /  >  sin(kLT/2) 

77      AT^  C°S  k(z"Zm)  dS      =     kLT/2  (1  -33) 

Measured  voltage  outputs  must  be  corrected  by  this  factor  before  the 
results  are  used  in  the  calibration  equation  when  transducer  dimensions 
are  an  appreciable  fraction  of  a  wavelength.  For  the  LC-10's  used  in 
our  experiments,  this  factor  is  0.05  dB  at  1kHz;  for  the  LC-32's,  it  is 
0.27  dB. 

3.      Limits   of   Derivation 

The  derivation  of  the  compliant  tube  reciprocity  parameter  has 
been  done  under  four  important  simplifying  assumptions.  To  the  degree 
that  these  approximations  do  not  represent  the  real  system,  the  results 
obtained   by   absolute  calibration    will    be   in    error. 

The  first  of  these  approximations  is  the  use  of  the  solutions 
to  the  lossless  wave  equations  in  computing  the  total  energy  in  equation 
(1.28).  The  PVC  wall  does  attenuate  the  sound  wave,  resulting  in  a 
resonance  sound  field  which  resembles  that  derived  by  Kinsler,  Frey, 
Coppens,  and  Sanders  [Ref.  10: pp .  206-209]  .  Figure  1.1  shows  the 
measured  variation  of  the  pressure  amplitude  with  depth  for  the  n=4 
mode.  The  pressure  amplitudes  at  separate  maxima  will  therefore  not 
be  identical,  as  is  required  by  our  development,  which  assumes  an  infi- 
nite  standing   wave   ratio. 

A  second  approximation  is  the  assumption  of  linear  behavior 
for  the  polyvinyl  chloride  wall.  In  fact,  PVC  has  no  elastic  region  for 
strain,  as  was  discussed  in  the  derivation  of  the  wall  motion  equation. 
The  nonlinearity  of  the  wall  permits  mode  coupling  into  other  resonance 
modes,  particularly  when  two  resonances  have  spatial  pressure  maxima 
near  the  same  position  and  are  related  subharmonically .  The  result  of 
such  coupling  when  the  column  is  driven  at  one  frequency  is  to  depress 
the  amplitude  of  individual  pressure  maxima  by  an  amount  proportional 
to  the  degree  of  coupling  at  that  position;  again  violating  the  require- 
ment  for   equal    pressure  amplitudes   at  the  two   test   positions. 
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Third,  the  reciprocity  concept  requires  that  the  receiving 
transducer  approximate  an  infinite  mechanical  impedance  to  the  medium. 
While  piezoceramics  do  meet  this  approximation  very  well,  the  soft 
mounting  material  used  the  the  hydrophone  body  and  gas  bubbles  on 
the  hydrophone  surface  do  not.  Their  presence  could  significantly 
distort  the  sound   field    in   the  vicinity  of  the  transducer. 

Finally,  the  projector  used  to  drive  the  column  into  resonance 
does  not  produce  a  sound  field  identical  to  that  generated  by  a  small 
piezoelectric  transducer  suspended  in  the  water  column.  The  local 
sound  field  near  the  cylinder  is  some  combination  of  radial  modes  which 
match  the  cylindrical  boundary  conditions  to  the  overall  plane  wave, 
coupling  the  driving  transducer  to  the  standing  wave  system.  The 
effects  of  this  coupling  factor  are  slightly  frequency-dependent,  and 
are   not  accounted   for   in   this   model. 

Our  measurements  demonstrate  that  the  simple  model  derived, 
despite  these  deficiencies,  produces  consistent  results  to  within  ±0.8 
dB,  but  does  not  explain  differences  from  a  standard  comparison  cali- 
bration of  up  to  2.1  dB.  The  root-mean-square  average  difference  is 
1.5±0.4  dB  for  the  range  of  hydrophone  types  and  frequencies  tested. 
A  method  and  apparatus  for  performing  this  type  of  calibration  experi- 
ment,   and   the   results   of  our  experiments,    follow. 


21 


Figure   1.1         Pressure   Amplitude   Depth    Profile    (n=4) 
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II.    COMPLIANT  TUBE   RECIPROCITY   EXPERIMENTS 

A.      RESONANCE   PARAMETER  MEASUREMENT 

The  formula  for  the  reciprocity  parameter  derived  in  the  preceding 
section  requires  the  accurate  measurement  of  three  resonance  parame- 
ters: the  resonance  frequency,  the  sound  speed  at  that  frequency, 
and  the  resonance  quality  factor:  Each  of  these  measurements  involves 
finding  an  extremum  from  a  set  of  .data  samples,  which  will  inevitably 
contain  noise.  The  experimenter  who  attempts  to  resolve  extrema  from  a 
discrete  sample  set  by  comparing  measured  signal  outputs  performs  the 
discrete  equivalent  of  differentiation,  which  enhances  noise  and 
increases    net  error   in   the   result. 

The  solution  to  accurate,  reproducible  resonance  parameter  meas- 
urement which  D.  Conte  and  S.  Garrett  originally  developed  [Ref.  12] 
combines  the  high  data  rate  of  automated  measurements  with  the  power 
of  curve  fitting.  The  intuitive  explanation  of  the  effectiveness  of  this 
method  is  the  human  analogue;  rather  than  take  discrete  measurements, 
experimenters  instead  "swing"  the  system  slowly  through  the 
extremum,  and  let  the  eye  interpolate  the  location  of  that  point.  The 
optimum  automated  processing  analogue  for  this  method  is  curve  fitting 
to  the  minimum  expansion  of  any  function  y(x)  which  will  resolve 
extrema,    the   second   order   equation 


y(x)   =   C0   +   C,X    +   C2X2  r2     -|) 


The  three  coefficients  of  this  curve  can  be  found  for  any  set  of  three 
or  more  measurements;  if  more  than  three  measurements  are  made,  the 
nonlinear  regression  methods  described  by  Bevington  [Ref.  13]  can  be 
used  to  reduce  the  noise  error  in  the  measurement  of  the  extremum,  as 
long   as   the   span   of  the   sample    set   around   the   extremum    is    small    enough 
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that  equation  (2.1)  is  a  valid  approximation.  The  regression  equations 
for  the  three  coefficients  are  algebraic,  and  are  derived  in  Appendix  B 
(equations  (B.12),  (B.13),  and  (B.14)).  The  location  of  the  extremum 
(again  from  Appendix  B)  is  displaced  from  the  coordinate  origin  of  the 
sample  set  by 


-  x0  =  —  Ci/2Cj 


(2.2) 


In  pur  application  of  this  method  to  measurement  of  the  resonance  peak 
frequency  ,  the  sample  set  is  composed  of  seven  hydrophone  output 
voltage  amplitudes  measured  at  equal  frequency  intervals;  the  sample 
set  origin  is  shifted  until  the  estimated  peak  frequency  falls  within  one 
interval  of  the  center  frequency  of  the  sample  set.  This  is  the  peak 
frequency  locating  algorithm.  Figure  2.1  is  a  example  of  this  technique; 
the  asterisks  mark  data  points  used  to  establish  the  coefficients  for  the 
curve  and  to  compute  the  location  of  the  maximum,  and  points  marked 
with  circles  are  discarded  from  the  final  fit  because  they  are  more  than 
three  steps  from  the  extremum.  The  sound  speed  is  found  by 
sampling  the  pressure  amplitude  at  regular  intervals  in  depth  from  the 
surface,  and  again  shifting  the  set  used  for  curve  fit  until  the  esti- 
mated maximum  or  minimum  is  within  one  step  of  the  center  of  the  set. 
The  positions  thus  located  occur  at  quarter  wavelength  steps,  and  can 
be  fit  to 


zm  =  mX/4  +  z0  (2.3) 

to  obtain  the  wavelength;  the  sound  speed  is  then  just  the  product  of 
the  resonance  peak  frequency  and  the  wavelength.  This  is  the  sound 
speed   measurement   algorithm. 

The    curve    fit    in    frequency    also    yields    a     reduced-noise    value    for 
the      peak     amplitude     at      resonance,      which      is        essential     to     accurate 
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measurement  of  the  half-power  frequencies  required  for  the  Q  measure- 
ment algorithm  (equation  (1.24)).  Once  the  half-power  amplitude  is 
known,  the  half-power  frequencies  can  be  found  by  sampling  outward 
from  the  peak  frequency  until  a  pair  of  samples  brackets  this  value  at 
the  minimum  frequency  interval  permitted  by  the  equipment  (in  our 
case,  1  Hz).  The  half  power  frequency  is  then  calculated  by  linear 
interpolation. 

The  curve-fitting  method  can  resolve  extrema  in  noise  to  a  degree 
far  surpassing  the  eye  integration  method  it  was  drawn  from.  It  has 
the  additional  advantages  of  sound  mathematical  basis,  speed,  reprodu- 
cibility, and  independence  from  observer  biases.  It  also  has  the  consid- 
erable advantage  of  being  easier  on  the  observer.  The  method's  main 
drawback  is  its  bias  toward  finding  the  parabolic  curve  in  any  data, 
including  pure  noise;  care  must  be  taken  in  selecting  stepsizes  to 
ensure  that  the  changes  in  the  signal  are  not  so  small  that  they  are 
overwhelmed  by  noise,  and  to  ensure  that  the  sample  set  span  is  not 
larger  than  the  region  about  the  extremum  where  a  second  order 
approximation  is  valid.  The  best  method  for  checking  new  cases  is 
graphical  display  of  the  data  set  and  the  curve  for  the  operator 
approval . 

B.      DESCRIPTION   OF   EQUIPMENT 

1.      Control    System    (Figure  2.3) 

The  experiment  is  controlled  by  a  venerable  Hewlett  Packard 
HP-85  laboratory  computer  with  32K  memory  extension,  IEEE-488  bus 
controller,  and  appropriate  ROM  extensions.  The  program  used  is 
documented  in  Appendix  C.  Commands  generated  by  the  program  over 
the  data  bus  control  the  experiment's  configuration  changes  during 
measurements  through  the  two  dual  remote  coaxial  switches  (Hewlett 
Packard  HP  59307)  and  a  Hewlett  Packard  HP  3421  Data  Acquisition 
Control  Unit.  The  latter  unit  is  configured  with  two  actuator  relays 
and  eight  measurement  channels  (option  020).  The  two  actuator  relays 
control   the  operation   of  the   reversing    DC   manipulator   drive  motor. 
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Figure   2.1         Parabolic   Curve    Fit   in    Frequency 
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Figure  2.2       ACTRE   Test   System 
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Figure  2.3        Control    System    Schematic 
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The  manipulator's  screw  drive,  which  is  described  below,  also 
drives  a  Hewlett  Packard  HED  680  position  encoder  that  increments  a 
Hewlett  Packard  HP  5016  universal  counter  in  totalizing  mode.  The 
changes  in  this  totalized  count  provide  (via  the  data  bus)  the  program 
with   a   measurement  of  actual   manipulator  motion. 

2.  Transmitter   (Figure  2.4) 

A  Hewlett  Packard  HP  3314  function  generator  provides  a 
stable  bus-controlled  sine  input  signal  at  two  voltage/current  levels, 
and  a  synchronization  signal  at  the  selected  frequency  and  constant 
voltage.  The  high  voltage  output  is  amplified  and  used  to  drive  the 
piezoceramic  transducer  at  rms  voltages  of  10  to  20  volts  AC.  The 
low-voltage,  high-current  demands  of  the  USRD  J11  projector  [Ref.  14] 
are  met  by  low-gain  power  amplification  of  the  function  generator's  low 
voltage  output.  The  transmit  remote  coaxial  switch  (714) 2  selects  the 
signal  source  from  these  two  options  and  feeds  the  signal  to  the 
reversible  transducer  switch  (706A),  to  the  J11  projector,  or  through  a 
test  attenuator  which  returns  a  known  signal  to  the  receiver  input 
switch    (706B). 

A  Hewlett  Packard  HP3478  digital  multimeter,  configured  as  a 
remotely  controlled  5  1/2  digit  AG  voltmeter,  is  connected  to  the 
transmit   signal    line  to  measure  the  transmitting    voltage. 

3.  Receiver    ( Figure  2.4) 

The  test  transducers  used  in  these  experiments  are:  two 
Celesco  LC-10  piezoceramic  hydrophones  (Figure  2.5),  selected  for  their 
small  size  and  good  transmitting  response;  and  a  Celesco  LC-32  piezo- 
ceramic hydrophone  (Figure  2.6)  of  larger  size.  The  voltage  sensitivi- 
ties of  all  three  units  were  measured  by  comparison  with  a  reference 
standard  USRD  H56  hydrophone  [Ref.  14].  The  results  are  listed  in 
Table  II.  The  LC-10's  are  mounted  on  0.375  inch  thinwall  copper 
tubing;    the    LC-32   is    mounted    on    0.4375    inch    tubing.    Signal    cables    are 


2The    remote     switches    are     identified     here     by    their     IEEE-488    bus 
address    codes. 
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TABLE   II 
Comparison   Calibration    Resul-ts 

(All   values   ±   0.8  decibels) 


Frequency 
(Hz) 


740 

860 

1000 


LC-10#2166 
&  2m  cable 

LC-10#2338 
&  9m   cable 

LC-32#85 
&   11m   cable 

-206.5 

-206.8 

-212.3 

-206.8 

-207.3 

-213.4 

-206.9 

-208.8 

-214.3 

routed  up  the  interior  of  the  tubing  to  minimize  effects  on  the  sound 
field.  The  two  test  transducers  are  suspended  from  the  manipulator 
platform.  Their  output  is  returned  from  the  top  of  the  manipulator 
assembly  by  two  of  three  coaxial  cables  of  the  same  length  and  meas- 
ured capacitance.  The  third  of  these  cables  connects  the  test  attenuator 
into  the  receiver  input  switch  (706B),  which  makes  all  receive  signal 
path  wiring  nearly  equivalent.  The  reversible  transducer  switch  (706A) 
selects  either  its  receive  or  transmit  lines;  the  receiver  input  switch 
(706B)  selects  the  lockin  analyzer  input  from  the  reversible  transducer 
receive   line,    the  test   hydrophone,    or  the  test   attenuator. 

The  EG&G  PAR  model  5204  lockin  analyzer  filters  out  the 
components  of  the  received  signal  which  are  in-phase  and  in-quadrature 
with  the  reference  (synchronization)  signal.'  Their  Pythagorean  sum  is 
used  to  obtain  the  magnitude.  The  analyzer  is  typically  operated  with 
an  integration  time  of  300  milliseconds,  which  corresponds  to  an  equiva- 
lent noise  bandwidth  of  0.8  Hz;  this  causes  the  lockin  analyzer  to  act 
as  a  very  sharp  filter  synchronized  to  the  transmitter  signal.  The 
analog  output  of  the  model  5204  is  passed  to  a  measurement  channel  of 
the   HP  3421    for  digitization. 

4.      Mechanical    System 

The  tests  are  conducted  in  a  6-inch  (154  mm)  inside  diameter 
Schedule  40  polyvinyl  chloride  tube  mated  to  the  active  surface  of  a 
USRD     J11     electrodynamic     underwater     projector     [Ref.    14].        The     tube 
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Figure  2.5        Celesco   LC-10   Hydrophone 
(Reprintecf  with    permission   of   Celesco    Industries,    Inc.) 


stands  1180  mm  above  the  Jll's  active  surface,  and  is  filled  to  a  depth 
of  about  1  meter;  its  wall,  which  is  smooth  on  the  interior  except  for 
small  regular  extrusion  marks,  is  8  mm  thick.  The  configuration 
resembles  the  USRD  G40  calibrator  used  for  comparison  calibrations 
[Ref.  14].  (This  system  was  originally  constructed  by  Mills  and 
Garrett  for  their  work  [Ref.  15]  with  fiber  optic  gradient  hydro- 
phones.) Equalization  for  the  J11  driver  is  provided  by  air  pressuriza- 
tion  of  the  installed  air  bag,  monitored  with  a  water-filled  U-tube 
manometer. 

The  test  transducers  are  suspended  on  thinwall  tubing 
supported  on  a  rigid  table  which  is  raised  and  lowered  on  two  one-inch 
rod  tracks  by  a  motor-driven  screw  drive;  maximum  range  of  motion  is 
800  mm.  The  drive  controls  were  described  in  the  paragraphs  on  the 
control    system   above. 
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Figure  2.6        Celesco    LC-32   Hydrophone 
( Reprinted  with    permission   of   Celesco    Industries,    Inc.) 


C.      PROCEDURES    FOR   EXPERIMENT 

1  .       Preparation 

Before  the  test  transducers  are.  suspended  in  the  water 
column,  the  experimenter  must  measure  the  reversible  transducer's 
input  capacitance,  so  that  the  transducer  current  can  be  computed  from 
the  transmitter  input  voltage  and  drive  frequency.  The  active  surface 
lengths  of  both  hydrophones  should  also  be  measured  for  use  in 
computing  the  transducer  amplitude  averaging  factor  (equation  (1.33)). 
The  test  transducers  are  then  suspended  from  the  manipulator,  wetted 
with  a  detergent  solution,  immersed  in  the  water  column,  and  brushed 
to  remove  bubbles  from  the  active  surface.  If  necessary,  the  tube's 
inner  wall  should  also  be  brushed  free  of  bubbles.  The  water  column 
depth    is   measured   with    a   meter   stick. 
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2.      Mode  Selection 

Two  parameters  are  available  to  vary  the  resonant  frequency 
of  the  water  column  at  which  the  calibration  is  done:  the  water  depth 
L  and  the  mode  number  n  (equation  (2.4)).  The  procedure  requires 
wavelengths  much  larger  than  the  coaxial  lengths  of  the  transducers. 
In  particular,  the  wavelength  should  be  much  longer  than  the  coaxial 
length  of  the  reversible  transducer  T  to  ensure  that  T's  surface  is 
entirely  within  the  velocity  node  and  feels  little  or  no  impressed  pres- 
sure when  T  is  transmitting.  The  wavelength  should  also  be  at  least 
four  times  the  length  of  the  test  hydrophone  H.  Practical  experience 
has  shown  that  calibration  results  are  most  accurate  when  T  and  H  are 
positioned  at  two  different  pressure  maxima  of  the  n=3  or  higher  order 
resonance   modes. 


L  =  -(n  +  1/2)=-  (n  +  1/2)  (2.4) 

2  2f 


3.       Resonance   Parameter  Measurements 

After  positioning  the  test  transducers  at  the  approximate  loca- 
tions of  the  deepest  and  next  deepest  accessible  spatial  pressure  maxima 
with  similar  pressure  amplitudes,  the  operator  uses  the  peak  frequency 
measurement  algorithm  to  find  the  resonance  frequency  f.  The  trans- 
ducers   are  then    withdrawn    above  the  first   pressure  minimum. 

One  of  the  transducers  is  now  used  to  measure  the  locations 
of  maxima  and  minima  in  the  standing  wave  pattern.  Linear  regression 
of  this  data  gives  the  quarter  wavelength,  which  in  turn  gives  the 
sound    speed. 

Finally,  the  transducers  are  precisely  repositioned  at  the 
maxima  where  the  resonance  frequency  measurement  was  made.  A 
precise  measurement  of  the  peak  frequency  and  of  the  Q  factor  are 
made,  completing  the  set  of  measurements  needed  to  compute  the  reci- 
procity   parameter  from   equation    (1.32). 
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•     4.      Calibration   Measurements 

Before  making  the  final  measurements,  the  experimenter 
switches  the  attenuated  test  signal  into  the  receiver  output  to  measure 
the  net  effects  of  cable  losses,  preamplifier  gain  drift,  and  deviations 
between  the  lockin  analyzer's  and  the  digital  voltmeter's  measurements. 
This  correction  is  applied  along  with  the  amplitude  averaging  correction 
factors  from  equation  (1.33)  to  the  calibration  voltage  measurements 
that  follow. 

The  transmitter  is  switched  to  the  reversible  transducer  T, 
and  the  receiver  to  the  test  hydrophone  H.  The  transducer-to- 
hydrophone  output  voltage  E(TH)  is  measured,  corrected,  and 
recorded,  as  is  the  transmit  voltage  V(T)  .  The  latter  is  used  with  the 
capacitive  admittance  of  T  to  compute  the  transmitter  current    l(T). 

The  transmitter  is  then  switched  to  the  projector,  and  the 
drive  voltage  adjusted  to  give  a  sound  field  of  the  same  amplitude  at  H 
as  was  used  in  the  previous  step.  The  projector-to-hydrophone  output 
voltage  is  recorded  (it  is  equal  to  E(TH),  within  the  limits  of  the 
projector  drive  adjustment  error).  The  final  calibration  measurement  is 
made  with  the  receiver  switched  to  T  with  the  projector  drive 
unchanged;  the  result  is  the  projector-to-transducer  output  voltage 
E(PT). 

These  measurements  complete  the  requirements  to  compute  the 
test  hydrophone  voltage  sensitivity  from  equation  (1.5).  The  same  set 
of  measurements  can  be  used  to  compute  the  reversible  transducer's 
sensitivity   using 


MT_JEjr  (2-5) 


The    results    for    M    are    in    volts    per    pascal;    the    customary    units    quoted 
are  decibels    re   1V/uPa,    which    is   found   by   calculating 


Jt  =  20  log  M  -  120  p   6) 
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5.      Self-consistency   Check 

When  both  transducer  can  be  used  as  reversible  transducers, 
the  results  of  the  reciprocity  calibration  can  be  checked  promptly  by 
exchanging  the  transducers'  roles  in  the  calibration  measurements. 
Repeat  measurement  of  the  resonance  parameters  and  transducer  reposi- 
tioning   is    not   required. 

Unless  the  two  transducers  have  identical  nonlinear  responses, 
this  second  set  of  readings  will  confirm  the  accuracy  of  the  assumption 
of  reciprocity.  It  will  not,  however,  measure  the  error  in  computing 
the,  reciprocity   parameter  J,    which    must   be  estimated    independently. 

D.      RESULTS   OF   EXPERIMENTS 

1  .      Reciprocity   Tests 

Tests  were  run  on  the  three  transducers  in  all  pair  combina- 
tions at  the  three  accessible  modes:  n  =  3,  4,  and  5.  Modes  below  n=3 
were  inaccessible  because  the  LC-10  transducers'  transmitting  response 
was  too  small  for  reliable  measurements;  modes  above  n=5  were  attenu- 
ated into  the  'tails'  of  the  lower  resonances,  which  obscured  the  half- 
power  points  needed  to  measure  Q.  Comparison  of  a  number  of 
reciprocity  calibration  test  results  showed  reproducible  results  to  within 
0.8  decibels.  However,  the  calibration  results  in  Table  IV  disagree 
systematically  with  the  values  obtained  for  the  same  hydrophones  by 
comparison  with  the  USRD  standard  hydrophone  (Table  II)  by  1.5  deci- 
bels   (Table  V). 

2.      Error   Analysis 

The  self-consistency  error  quoted  in  Table  IV  was  determined 
by  reversing  the  electrical  roles  of  the  two  transducers  without 
changing  their  positions;  the  difference  in  computed  results  for  the  two 
measurements  is  the  degree  to  which  the  system  is  not  reciprocal.  The 
estimated  error  is  the  combination  of  the  self-consistency  error  with  the 
error   in    measuring   the    reciprocity   parameter,    listed    in    Table    III. 
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The  primary  source  of  self-consistency  error  is  the  difference 
in  pressure  amplitudes  at  the  two  maximum  locations  used,  which  we 
believe  is  caused  by  attenuation.  By  careful  selection  of  modes  and 
locations,  we  were  able  to  reduce  this  error  to  0.2  dB.  Another  lesser 
source  of  self-consistency  error  is  the  changing  effect  of  the  support 
tubing  used  for  the  transducers  on  the  sound  field  as  they  are  lowered 
into  the  water  column.  Insertion  of  the  transducers  and  their  supports 
into  the  water  column  displaces  about  0.4  percent  of  the  water  column 
in  volume,  causing  a  similar  change  in  the  column  depth  and  in  the 
positions  of  maxima  and  minima.  The  change  in  the  minima  is  the  most 
significant  error  encountered  in  measuring  the  wavelength  of  the 
standing  wave;  because  the  cosine  function  is  flat  to  first  order, 
however,  the  effect  on  transducer  positioning  at  maxima  is  negligibly 
small . 

Compliant  materials  associated  with  the  signal  cables  and  the 
transducers  appear  to  degrade  the  reciprocal  nature  of  the  geometry. 
This  error  is  reduced  to  very  low  values  by  selection  of  transducer 
location  and  by  eliminating  compliant  materials  (signal  cables,  for 
example)    in    contact  with   the  water. 

Parameter  measurement  error  is  dominated  by  the  apparent 
variation  in  the  system  quality  factor  over  the  length  of  the  tube.  An 
increase  of  6-10  percent  in  the  value  of  Q  determined  by  the  half-power 
frequency  method  was  observed  as  the  depth  of  the  measurement  was 
increased  by  one  wavelength;  our  model  for  the  compliant  tube  system 
does  not  account  for  such  effects,  and  is  therefore  limited  in  the  accu- 
racy with  which  it  predicts  the  reciprocity  parameter  to  about  0.8  dB. 
This  is  the  largest  recognized  source  of  measurement  error  in  the 
compliant  tube    reciprocity   method. 

The  systematically  lower  sensitivities  measured  by  the  reci- 
procity method  relative  to  those  measured  by  comparison  with  the  H56 
standard  are  not  yet  understood.  This  discrepancy  may  rise  from  the 
inadequacies  of  our  model  for  the  system  energy  balance,  from  the 
method  used  to  measure  Q,  or  even  from  technical  difficulties  with  the 
comparison    calibration. 
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3.      Measurement  of  Wall   Motion 

Since  sound  speed  can  be  measured  directly,  it  is  not  neces- 
sary to  attempt  the  difficult  measurement  of  the  actual  elastic  strain  in 
the  tube  wall  as  part  of  the  calibration  process.  However,  the  degree 
of  correspondence  of  the  actual  strain  behavior  to  the  model  derived  in 
the  first  chapter  is  an  important  indication  of  the  validity  of  our 
method . 

The  circumferential  wall  strain  was  measured  by  attaching 
constantan  .250  inch  strain  gages  [Ref.  16]  to  the  exterior  surface  of 
the,  tube  at  the  depths  of  some  of  the  pressure  maxima  located  by  the 
test  transducer.  A  second  lockin  analyzer  used  as  the  sensing  element 
of  a  resistive  DC  bridge  measured  the  magnitude  and  phase  of  the 
strain  oscillations  sensed  by  the  strain  gage;  by  comparing  the  output 
of  the  strain  gage  with  the  magnitude  and  phase  of  the  pressure  meas- 
ured by  the  transducer  at  that  depth,  we  computed  the  complex  impe- 
dance (W,9)  of  the  wall  ring  (equation  (1  .21 ))  .  The  results  in  Table  VI 
agree  in  general  with  the  measured  sound  speed  in  the  tube.  Also,  the 
indicated  damping  phase  effects  are  reasonably  small,  justifying  our 
approximation  for  the  wall  potential  energy  computation  of  equation 
(1.23). 
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III.    UTILITY   OF   COMPLIANT  TUBE   CALIBRATION 

A.      LIMITS  OF  COMPLIANT  TUBE  CALIBRATION 

The  compliant  tube  reciprocity  method  can  be  extended  to  low 
frequencies  by  simply  increasing  the  length  of  the  water  column 
enclosed  in  the  PVC  tube,  or  by  decreasing  the  wall  thickness  to  diam- 
eter ratio  (equations  (1.11)  and  (1.20)).  The  practical  lower  frequency 
limi.t  is  the  length  necessary  to  support  the  n=3  standing  wave,  which 
from   equation    (2.4)    is 


f      -Z-£ 


The  upper  frequency  limit  depends  both  on  the  onset  of  azimuthal 
modes,  and  on  the  increasing  attenuation  of  the  sound  wave  by  the  PVC 
wall.  The  former  constraint  breaks  the  plane  wave  symmetry  when  the 
frequency  is  greater  than  (0.586#L/R)  times  the  fundamental  longitu- 
dinal mode  [Ref.  17].  The  attenuation  problem  causes  excessive  varia- 
tion of  the  measureable  system  quality  factor  with  depth  and  reduces  Q 
to  the  point  where  half  power  frequencies  are  obscured  by  adjacent 
resonances . 

Test  transducer  geometry  is  restricted  to  shapes  that  are  small 
enough  in  cross  section  that  reflections  from  their  lower  surface  do  not 
significantly  affect  the  standing  wave.  The  length  of  the  reversible 
transducer,  as  discussed  in  the  procedure  description,  should  be  as 
short  as  possible  consistent  with  adequate  transmitting  response.  (A 
flat  disk  transducer  might  be  more  effective  than  the  cylindrical  LC-10 
we  used.)  The  test  hydrophone  can  be  much  larger,  if  the  transducer 
amplitude   averaging   correction    is    properly   computed. 

Accuracy  of  the  method  is  primarily  limited  by  the  effects  of  the 
wall    attenuation.       A    derivation    which     accounted    for    these    effects,    and 
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for  the  lesser  effects  of  disturbances  to  the  standing  wave  caused  by 
the  test  transducers  and  their  support  tubing,  could  conceivably 
achieve  much    lower  calibration    uncertainties . 

A  significant  advantage  of  the  compliant  tube  method  is  the  use  of 
pressure  antinodes  as  virtual  rigid  boundaries  at  which  to  locate  the 
test  transducers.  Many  of  the  difficulties  associated  with  mounting  a 
transducer  in  an  real  boundary  which  must  have  some  compliance  are 
•thus   avoided. 

B.      POTENTIAL  APPLICATIONS 

Two  applications  for  the  compliant  tube  system  on  a  practical  scale 
have  been  suggested  to  the  author:  a  field/small  lab  calibrator,  and  a 
large,    very    low  frequency   calibrator. 

The  portability  and  flexibility  of  the  test  system  used  for 
compliant  reciprocity  recommends  the  method  for  general  calibration  use 
in  small  laboratories  and  in  the  field.  The  1 -meter  PVC  tube  provides 
calibration  capability  in  the  1  kHz  range;  a  2-meter,  10-inch  diameter 
tube  would  be  more  versatile  and  could  reach  frequencies  of  500  Hz. 
The  manipulator  assembly  need  not  be  as  complicated  or  massive  as  the 
one  used  in  our  experiment;  any  system  which  can  position  the  trans- 
ducers within  a  few  millimeters  is  adequate.  Because  the  method 
involves  resonance  measurement  methods  as'  well  as  the  techniques  of 
reciprocity  calibration,  it  may  also  be  of  interest  for  advanced  instruc- 
tional   laboratories. 

On  a  larger  scale,  a  10-meter,  30-inch  diameter  tube,  which  would 
fit  upright  in  a  two-story  open  bay,  could  reach  frequencies  below  100 
Hz  and  calibrate  units  similar  in  size  to  the  Navy's  TR-155  transducer. 
If  current  interest  in  very  low  frequency  underwater  acoustics  indicates 
a  need  for  an  alternative  to  the  traveling-wave  tube  system  described 
in  [Ref.  8],  compliant  tube  reciprocity  may  provide  the  answer  at  a 
very    low   cost. 

There  is  also  no  reason  we  are  aware  of  that  the  benefits  of 
compliant    waveguide    usage     should     be     limited    to     resonant     calibration. 
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Construction  of  a  compliant  traveling  wave  tube  calibrator  with  a  non- 
reflecting  termination  or  of  a  Helmholtz  resonant  calibrator  using  two 
compliant  tubes  should  also  be  practical.  The  Helmholtz  calibrator 
concept  in  particular  has  the  potential  of  extending  the  range  of  accu- 
rate absolute  calibration  for  the  underwater  acoustician  to  extremely  low 
frequencies. 
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APPENDIX  A 
GLOSSARY 

Subscripted  variables  in  the  text  are  represented  using  paren- 
theses,   due  to  word   processor   limitations. 

A  Cross   sectional   area  of  the  water  column   at   rest. 

A(T)  Active   surface  area   of  the  transducer. 

a  Instantaneous   cross   sectional   area   of  the  water  column. 

C(0  Parabolic  curve  fit   coefficients. 

c  Sound    speed    in    a   compliant  tube. 

D  Distensibility  of  the  boundary   wall. 

E(TH)  Transducer-driven,    test   hydrophone  output   voltage. 

E(PH)  Projector-driven,    test   hydrophone  output   voltage. 

E(PT)  Projector-driven,    reversible  transducer  output   voltage. 

f  Resonance  frequency. 

H  The  test  hydrophone  used  only  to  receive  during  the  reci- 
procity  calibration. 

h  Thickness   of  the  tube  wall. 

I  (T)  Transmitter  current   in    reversible  transducer  T. 

J  Reciprocity  parameter,  or  acoustic  transfer  admittance, 
between   two   points. 

K  Adiabatic   compressibility   of  water. 

k  Wave   number   at   resonance  =   2ir/\. 

L  Depth   of  the  water  column. 

L(T)  Transducer's   active   length    normal    to   the   plane  wave. 

(VI  Hydrophone   voltage   sensitivity. 

m  Number  of   quarter  wavelengths    below   the  water   surface. 

n  Resonance  mode   index. 

n  Unit   vector   normal    to   surface  element  ds. 

p  Acoustic  wave  excess   pressure. 

P  Pressure  amplitude  of  the  acoustic  wave. 

Q  Resonance  quality   factor. 

R  Inner    radius    of   the   tube. 
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R(m)  Damping   coefficient  for  tube  wall   vibrations. 

S  Transmitting   current   response. 

s  Surface  area   element. 

T  The    test    transducer    used     both    as     a     projector    and     as     a 

hydrophone  during  the   reciprocity   calibration. 

U  Volume  velocity  of  the  acoustic  wave, 

u  Particle  velocity  of  the  acoustic  wave. 

V  Volume  element. 

W  Complex   mechanical    impedance  magnitude  of  the  tube  wall. 

x(0)  Center  coordinate  of  the  data   set   used   for  curve  fitting. 

x(m)  Coordinate  of  extremum  of  the  curve. 

Y  Elastic  coefficient  for  tube  wall   expansion. 
y(i)  Sample  set  to  be  curve  fit. 

Z(E)  Transducer  blocked   electrical    impedance. 

Z(m)  Transducer  open-circuit   mechanical    impedance. 

z  Depth   from   surface  of  water  column. 

z(0)  Depth   measurement  offset. 

z(m)  Depth    of  the  mth    pressure  maximum/minimum. 

5  Sampling    interval   for  curve  fit   sample   set   y(i). 

s  Areal    strain    in   the  wall.      See   equation    (1.19). 

X  Length   averaging   correction   factor. 

8  Complex   mechanical    impedance  phase   angle  of  the  tube  wall. 

X  Wavelength   of  the   standing   wave. 

P  Instantaneous   local   water  density. 

p  Ambient  water  density. 

$  Transformation   factor. 

X2  Sum   of   squared   errors    in    a   curve  fit   to   data. 

co  Angular  frequency   at    resonance. 

£  Areal    expansion   coordinate. 
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APPENDIX   B 
PARABOLIC   SMOOTHING   ALGORITHM 

The  algebra   of    deriving   formulas   for  the    best  choice  of  coefficients    of  a 
nonlinear  equation,    such   that  the   sum  of   squared   errors   described   by 


x2  =  Exf=  £ 


y(x<)  -  yi 


(B.i) 


is  minimized  ,  is  described  by  Bevington  in  [Ref.  13].  The  minimum 
for  the  function  defined  by  equation  (B.I)  is  defined  by  the  selection 
of  coefficients   which    solve 


§£  =0 

dC 


(B.2) 


where   C    is      any   of  the   coefficients   of  the   curve  to   be   fit     to  the   data. 

The      simplest    function    that    will    predict    maxima    and    minima     is    the 
quadratic   equation 


y(x)  =  C0  +  CiX  +  C2x: 


(B.3) 


The   single      extremum   of   this   curve  occurs   at 


xm  -  x0  =  -C,/2C3 


(B.4) 


with    respect    to    the    coordinate    origin    located    at    x.       Let    the    sample    set 
be   seven    measurements   of   the  function    y    sampled   at    regular   intervals    5 
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y.  =  y(x0  +  i<5)    ;    i=  -3,  -2,  -1,  0,  +1,  +2,  +3.  (B.5) 

Evaluating  the  best  fit  function's  partial  derivatives  with  respect  to  the 
three  parameters  for  such  a  sampled  data  set  gives  equations  (B.6) 
thru  (B.11);  summing  and  cancelling  yield  equations  (B.12),  (B.13), 
and  (B.14).  These  analytic  results  can  then  be  used  with  equations 
(B.3)  and  (B.4)  to  obtain  the  location  of  the  maximum  and  its  ampli- 
tude,   or  the   smoothed   value  for  the  amplitude  at  any   location. 

(B.6) 
X-s  =  9<52C2-  3<5C,+  Co-  y-3 

X-2  =  4<52C2-  2c5C,+  Co-  y-2  (B-7) 

X-]  =   <5  C2  —   <5Ci+  Co—  y-i 

X  o  =   Co- y  o 

X+i  =   62C2  +   <5C,+  C0-  y  +  i 

x+2  =  4(52C2+  2(5Ct+  C0-  y+2 

X+i  =  9c52C2  +  3<5C,+  Co-  y  +  3  (B.11) 


(B.8) 

(B.9) 

(B.10) 


(B.12) 


84(5  C2  =  5(y  +  3+  y_3)  -  3(y  +  ,+  y_,-  4y0  ) 

28(5C,  =  3(y  +  5-  y_5)  4  2(y+2-  y_2)  +  (y  +  I-  y.()  (B"13) 

21C0  =    7y0   +  6(y  +  t+  y_,)  +  3(y+2+y_2)  -  2(y+3+  y_3)  (B.14) 
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APPENDIX   C 
ACTRE:    HP-85   BASIC   PROGRAM   FOR    RECIPROCITY    EXPERIMENTS 

A.  WHAT  THE   PROGRAM   DOES 

The  Automated,  Compliant  Tube  Reciprocity  Experiment  (ACTRE) 
program  provides  the  operator  with  a  menu  of  soft  key  controlled 
"tools"  for  resonance  parameter  and  calibration  measurements.  Its 
desjgn  was  based  on  the  need  for  maximum  flexibility,  rather  than 
maximum    unattended   operation. 

All  of  the  measurement  "tools"  are  subroutines  called  from  the 
HP-85  function  keys  by  the  operator.  The  major  procedures  are:  the 
peak  frequency  location  algorithm  (KEY  1);  the  spatial  optimum  location 
algorithm,  which  also  computes  wavelength  and  sound  speed  (KEY  2); 
the  Q  measurement  algorithm  (KEY  5);  and  the  reciprocity  calibration 
measurement  sequence  (KEY  7).  Supporting  procedures  provide  the 
actual  control  statements  for  remote  switching,  function  generator 
setup,  manipulator  positioning  of  the  hydrophones,  receiver  output 
measurement,  and  graphics  output.  Any  procedure  may  be  terminated 
by  calling  the  process  clear  sequence  (KEY  8).  Comments  on  the 
detailed    structure  of  the   program   are   included   with    the   source   code. 

The    program   was    developed   to    fit    the   demands    of  the    experimental 
apparatus, the   experimenter,    and   the    research    goals    of  this   work   as   well' 
as    the    limits    of    the    host    computer    permitted.       Time    and    core    memory 
constraints     severely     limited    efforts     to     make     the    program     friendly     or 
accessible  to   the   uninformed. 

B.  TEST    EQUIPMENT    BUS    CONNECTIONS 

The  code  is  designed  to  use  the  IEEE-488  bus-controlled  test 
equipment  described  below.  All  of  the  primary  measurements  are 
made  with  bus-controlled  equipment,  except  for  the  hydrophone  output 
voltage  measurement  by  the  EGo-G  5204  lockin  analyzer,  whose  output 
signal   must   be   converted    using   the   HP  3421    D/A      unit. 
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Bus   assignment   numbers    used    in   the  program: 

705  HP7470   hardcopy   plotter. 

706  HP59307A    remote  controlled   switch: 

A   switch:    reversible  transducer; 
B    switch:    receiver   input. 

707  HP  3314  function   generator. 

709  HP    3421     Data    Acquisition/Control     Unit,     with     Opt    020 

board    in    slot  0.    Channel    assignments: 

0:    Manipulator  motor  drive   up   actuator 
1  :    Manipulator  motor  drive  down    actuator 
5:    5204   lockin   analyzer  magnitude   signal. 
714  HP   59307 A    remote  controlled   switch    (transmit   side). 

720  HP   5316    Frequency   counter,    in    count  totalize   mode. 

723  HP   3478    digital    multimeter,    configured    as    an    AC    volt- 

meter. 


VARIABLE   ASSIGNMENTS 
1  .       Defined    Functions 

FNH(Xl)       Length         averaging        correction        function         (equation 

(1.33)). 
FNK(KO)       Cycles     K0     through     the     indices     1,2,3.     Used     in     plot 

scaling   algorithm. 
FNL(XO)       Returns   decibel    value  of  X0  =   20   log(XO). 
FNM(XO)       Returns     voltage     sensitivity     in     dB     re     IV/uPa     (equa- 

tion(2.6)). 
FNR(XO)       Returns      value     of     X0     rounded      to     three      significant 

places;     required    for   compatibility    with    HP3314    function 

generator   amplitude   command. 
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2.      Variables 

A(7)  Amplitude    sampling   array,    used    for   curve  fitting    algo- 

rithms. 

A2  Amplitude     ratio        of     hydrophone      output     voltage      to 

transmitter  input  voltage.  Measured  by  the  Amplitude 
Read   procedure. 

A4  Half   power  amplitude      in    Q   measurement   procedure. 

A5  Lockin      analyzer     output      signal     to      D/A      unit      (709) 

channel  5;  ranges  from  0-1  V  to  indicate  0-full  scale 
magnitude   reading. 

A7  Resonance   peak   amplitude,    computed   from   curve  fit. 

A8  Preceding      measurement       storage--in       Q      measurement 

procedure. 

B(3)  Plot   scaling    parameter  array 

B1-B9  Plot   parameters 


C1-C5 


C6 
C7 
CS 


Sums     used     in     wavelength     measurement     least-squares 
procedure    [Ref.    13]. 

CI:    Sum   of   N9. 

C2:    Sum   of   N9*N9.       ' 

C3:    Sum  of   L9. 

C4:    Sum  of    L9*N9. 

C5:    Sum   of    L9*L9. 
The  quantity   A   in   the   least-squares   fit. 
Resonance   sound    speed    (m/s). 
Error   estimate  for   C7    (m/s). 


D1 
D2 
D3 
D5 
D6 


Water  column   depth(mm) 

Water  column   cross    section    (mm2) 

Water  column    volume(m3) 

Length   of   T    (mm)  . 

Length   of   H    (mm)  . 
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El  E(PT):    Output  voltage  of  T   ,    P     driven    (V). 

E2  E(PH):    Output  voltage  of   H    ,    P     driven    (V). 

E3  E(TH):    Output  voltage  of  H    ,    T     driven    (V). 

All   frequencies   are   in    hertz(Hz)  . 

FO  Center  frequency  of   sample   set. 

Fl  First  peak  detection   /    Lower   1/2   power  point. 

F2  Second   peak  detection   /    Upper   1/2   power   point. 

F3  Peak  frequency  estimate. 

F4  Lower  search    boundary. 

F5  Upper  search   boundary. 

F6  Drive    frequency;     argument    passed    to    Amplitude    Read 

Setup   procedure. 

F7  Resonance  peak  frequency. 

F8  Error   estimate  for   F7. 

F9  Stepsize. 

G2  Receiver  preamplifier  gain. 

G5  Lockin    analyzer  output   expand    setting. 

G7  Calibration   attenuator  value.' 

G8  Drive      boost      factor     for      preamp      calibration      test      in 
Calibration    routine. 

HI  Length    averaging   correction    factor  for   T. 

H2  Length    averaging   correction   factor  for   H. 
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10  General    use 

13  Size  of  plot  files  X(),    Y(). 

14  Selection    index   for  Operator    Input  procedure. 

15  Number   of    time    constants    T5    to    delay    before    measure- 
ments . 

18  Select   index   for  Switching   procedure. 

19  Wavelength   measurements   count. 

J3  Reversible  transducer(T)    transmit  current    (A). 

J7  Reciprocity  factor      (V*A/Pa2). 

K0  Water   density        (kg/m3). 

K3  Plot   scaling   array    index. 

L9  Pressure  maximum/minimum   position        (m)  . 

Ml  Hydrophone   voltage   sensitivity     of   T      (V/Pa)  . 

M2  Hydrophone   voltage   sensitivity   of      H      (V/Pa). 

N0-N2  Position   counts-    used   to   compute  manipulator  motion. 

N9  Number   of   1/4   wavelengths    from   surface    corresponding 
to    L9. 

Q7  Measured   quality   factor  of   resonance. 

R0  Position       encoder         count       ratio       for       down        motion 

(counts/mm)  . 

R1  Manipulator  drive    rate      (counts/sec),    down   motion. 

R2  Manipulator  drive      rate    (mm/sec),    down    motion. 

R3  Allowed   projector  drive   adjust  error. 
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R4  Allowed   peak  frequency   detection    range. 

R5  Allowed   min/max   detection    range. 

R7  Q   measurement  coarse   stepsize  factor. 

R8  Manipulator  up/down   motion    rate   ratio. 

The  S   variables   are  all   control   flags/indices. 

50  General    use. 

51  Read     enable  flag. 

52  Manipulator  motion    stop   flag. 

53  Plot  option   enable  flag. 

54  Fast  display    refresh   flag. 

55  Drive  voltage  select   index: 

0:    use  J11    projector   drive   voltage  V4; 
1:    use   piezoceramic  drive   voltage  V3 . 

56  Automatic   drive   voltage   control    enable   flag. 

57  Resonance   peak   detection    index. 

S9  Plot   execute   flag   and   abscissa    select   index. 

All   times   are   in   milliseconds. 

T2  Manipulator   motion    time   adjust. 

T5  Lockin    analyzer  time  constant. 

T7  Additional    delay   time   before   measurement. 

T8  Size  of   additional    delay   T7   to   be   used   when    required 

U0  Curve   fit   amplitude   at   center  of   sample   set. 

Ul  Curve  fit   slope*stepsize. 

U2  Curve  fit   curvature     *(stepsize)  2  . 

U9  Stepsize   used    in    last   curve  fit    (    for   plot   procedure). 
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All   voltages   are  in   volts(V)  rms,    except  where   indicated. 

VI  Transmitter  input  voltage. 

V2  Hydrophone  output  voltage. 

V3  Piezoceramic  drive  voltage(Vp-p) . 

V4  J11    projector  drive  voltage(Vp-p)  . 

V5  Lockin      analyzer      fullscale      setting;       input      as      milli- 

volts,and   converted   to  volts. 

V6  Function   generator  drive  voltage    (Vp-p)  . 

V7  Transmitter   input  voltage  at   resonance. 

V8  Curve      fit      amplitude      at       F7,       returned      by      5-point 

Voltage  Measurement   procedure. 

W0-W3  General    use  working   variables. 

X(50)  Abscissa   plot   storage   array    (Hz  or  mm). 

Y(50)  Ordinate   plot   storage   array    (mV   rms). 

XO  Reference   x-coordinate  for  curve  fit   plot    (Hz   or  mm). 

All    variables    beginning   with    "Z"   are   in    millimeters(mm)  . 

Z(7)  Sampling    positions   for   X   Measurement   procedure. 

ZO  Manipulator  table  position. 

Z2  Estimated    position   of   pressure   minimum/maximum. 

Z4  Position    of  center   of   sample    data    set/    operator-ordered 

position   to  move  to. 

Z8  Measured    stepsize   in    X   Measurement   procedure. 

Z9  Ordered   stepsize. 


3.      Character   String   Variables  . 

AS(128)  Title;    used   to   identify   printout   and    plots. 

MS(20)  Code   identifier. 

SS(5)  Remote   switch    position    identifier. 

WS(32)  General    character   input. 
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D.  SOURCE  CODE 


1000 

1010 
1020 

1030 
10  40 
1050 

1060 

1070 

1090 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1130 

1190 

1200 

1210 

1220 

1230 

12  40 

1250 

1260 

1270 

1280 

1290 

1300 

1  3  1  0 

1320 

'3  30 

1340 

1350 

1360 

i  3  7  0 

1330 

1390 

1400 

1  4  1  0 

1420  ! 

1430 

1  4  4  0 


program  actre-7 

Automated,  Compliant  Tube  Reciprocity  E: 

M.B.  Johnson,  Dec  84 

PRINTER  IS  2  d  OPTION  BASE  IS  DEG 
DIM  A$C128],H*[5],W$[32],S*[5J 
A$="ACTRE-7  mod  8     "  @  CLEAR  @  DISP  A* 
!  Functions 


peri  ment 


DEF 
DEF 
DEF 
DEF 
DEF 
Wl  =  . 


XI) 
Kl) 
XI) 
XI) 


=  X  1  /S  IN ( RTD  C  X  1 ) ) 

=  Kl    MOD    3+1 

=  20*LGT(Xi) 

=  FNKXD-120 


FNH 

FNK( 

FNL( 

FNM( 

FNR(Xl)  !  rounds  to  3  significant  figures  for 

01 


MP' 


14 


!  begi  n 
IF  IP  C X 1 
IF  IP(X1 

1  end  1  o 
FNR=M1*I 
FN  END 
1  end:  f 
1  Variab 
REAL  A  (7 
REAL  D3, 
REAL  NO, 
SHORT  B( 
SHORT  RO 
SHORT  X( 
INTEGER 
INTEGER 
1  i  n  i  t  i  a 
S1,S2,S6 
S3,S4,S9 
D2=PI«1S 
B2=l  8  3 
R0=512.5 
R3=. 01  ! 
R4. 


loop 

XI  THEN  X1=X1*1000  @  W i  ="- W 1/1000  §  GOTO 

) >9  THEN  Xl-Xl/10  @  Ml =M 1*10  @  GOTO  114i 

op 

P (X 1*100+. 5) 


1 1  4  0 


1 


14  =  6  @ 

!  4  ■-■  3  S 

14  =  2  @ 

R7=120 
B  t 

15  =  5  e 
KO-998 


)  ,  B  (  2 

T 
n 


ct  i  o 

i  ni 

A2,A 

,N2, 

,Bi, 
1  ,R2 
>  ,  X  0 
,B9, 

,S1, 
val  u 


ns 
ti 

4, 
E3 
UO 
B2 
,R 
|Y 
10 
S2 
es 


al  i  zat i  on 

A5  ,  A7 , A8 , C 1 , C2 , C3 , C4 , C5 , C6 , C7 , C3 , 

,  FO  ,  F  i  ,  F2  ,  F3  ,  F7  ,  G2  ,  J  3  ,  J  4  ,  J  7  ,  KO  , . .9 

M1  IP  V  !  V 1    '/' 8  Y 3  7  (  7 )  7  0  7  ->  T\ 
,33,B4,B5,B6,F4,F5,F6,F3,F9rG5,37 
3,R4,R5,R7,R8,T2,T5,T7,T8,U9,V3,V 
(50) ,2? 

,13,14, I5,I8,I9,K3,K?,N3,N4,N9 
,S3,S4,S5,S6,S7,S9 


CI, 
,11 

::. 

,33 
4  ,  V 


Z3 


151 
289 

Rl 

ad 

I 

UB 
0  @ 
UB 
ou 

2  § 

0  ! 


/4 

.  6 
*6 

ju 

ps 

20 

G 

2  0 

ti 

B 

r 

/m 


39 

00 

st 

ak 

60 

DSU 

6  0 

tep 

(3) 

ead 


mm"2 

@  68  =  30 

@  R2-.6/R0  3  RB=i  5  T 2  =  1 5 3 

crro  r 
detect  error 


move  -  - ;  i  j  I  .: 


B  206  0  '  position  input 
1  depth, temperature 

ratio 
=2.5  !  plot  scaling  parameters 


delay  factors 
water  density 
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1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 

I  CCA 

1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1630 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1730 


QU 
0U 
0U 
0U 
0U 
0U 
QU 
14 


13 
CI 
P7 

Ml 
V4 

S6 
18 

GO 


CL 

DI 
IF 
IF 
DI 
KE 
FO 
Ip 
NE 
3  4 
GO 


begi 
TPUT 
TPUT 
TPUT 
TPUT 
TPUT 
TPUT 
TPUT 
=  1  @ 
end: 
bsgi 

,19, 
,C2, 

iFB, 

,M2= 

=  .00 

=  1  ! 

=  1  e 

SUB 
and : 
bsgi 
hold 
EAR 
3P  U 

S3 

36 
SP  3 
Y  LA 
R  17 

3  4 
XT  I 
=  0 
TO  1 


n :  eq 

709 

707 

723 

707 

709 

720 

709 

GOSU 

squi 

n:  fr 

N9,N4 

C3,C4 

C7,C3 

10000 

01  9 

ADC 

GOSU 

1800 

f  res 

n  :  p  r  o 

here 


ui  p me 

"RS" 
H  pp  ii 

M  L)  T  I' 

"API 
"FIR 
"FN1 

"CL3 
9  206 
pment 
esh  s 
=  0  ! 
,C3,C 
,Q7,V 

00  ! 

V  3  =  10 
on 
B  312 

1  set 
h  si  a 
c  s  s  s  i 

■for 

8590 


n  t  setup 
!  cl ear 
1  cl ear 
!  xmi  t 

OMV" 

0RA021N5 

2"  !  pos 

5" 

0  !  LI  A 
setup  — 

late  -for 

plot  HI 

6=0  !  wa 

6  =  0 

0  dB  re 
3  T7  =  0 


A/D  uni  t 
3314  -fen  genera 
vol tmeter 


tor 


LS5;T1"  !  a/d  receiver  output 
ition  counter 

parameters 


processing  ====== 

e    2r.i1/4    file    size: 
velength    sums 

lV/uPa 


w 1 / 4  count 


0  I  switch 
main  menu 
te  for  processing"  —  —-' 
ng  stages -  =  =  --• -  =  =  ■=•--"=  = 
operator  menu  selection 


;  Af,F7, 

"plot  on 

"ADC  on" 

Z0=";Z0 


SING 

THEN  DISP 
THEN  DISP 

*r 

BEL 

=1  TO  49999 

THEN  1760  !  fast 

7 

660 


F3,C7,C8,Q7,FNM( 
";EL3E  DI3P  " 
ELSE  DISP  "ADC 


Ml  )  ,FNM(M: 

plot 

« i  r  'i 

Of  T 


BTt 


refresh  on  return  from 


.  e  y  r  o  i 


i  n : 


end:  processing  stage333-33 
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1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1380 
1390 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1990 
1990 
2000 
2010 
2020 
2030 
20  4  0 


ON 
ON 
ON 
ON 
ON 
ON 
ON 
ON 


begin: 
KEY* 

KEY* 
KEY* 
KEY* 
KEY* 
KEY* 
KEY* 
KEY* 
S4=l  !  i 
RETURN 
e  n  d :  m  a 
begin: 
p  r  o  v  i 
OFF  TIME 
OUTPUT  7 
51,52=1 
1 3 , S3 , 59 
ON  KEY* 
ON  KEY* 
ON  KEY* 
PRINT  "0 
54-1  !  f 
GOTO  165 
1  end:  c 


s  =  =  =  MENUS  =  =  =  — «  =  =  =  ' 
main  m8nu=======-==: 

l,"Fn"  GOSUB  2650 
2,CHR*(11)  GOSUB  4 1 5< 
3, "move"  GOSUB  7990 
4,"opr"  GOSUB  2060 
5,"Q"  GOSUB  3380 
6,"a/pt"  GOSUB  8120 
7,"CAL"  GOSUB  5190 
9,  "clear"  GOTO  1930 
ast  refresh 


men 
lear 
s  ex 

i  e 

;  "0 
flag 

!  P 

"mai 
"  w  i  p 
"qui 
RATO 
t  r  o 
en 


slate  men 
it/rsfrcsh 


u 

pro 
*  2 

PUT 


ess 
709 


si  n 


OFF  TIMER 
PNO"  (3  OUT 

reset 
lot  dump 
n"  GOSUB  1 
e"  GOTO  15 
t"  GOTO  85 
R  CLEAR/RE 
fresh 
t s r  processing  stage 


g 

OPNl 


GOO 
iO 
140 

:3£T' 
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2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2270 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 


3  INPUT  I 


INPUT  Dl ,W1 


;  =  =  =  =  =  =  proceDURES  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  -  =  =  =  =  =  =  =  =  =  =  =  — 

1  begin:  operator  input  procedure  (HP-HUMAN  bus) -  =  "  =  =-" : 

IF  14=0  THEN  2030  ELSE  2130 

!  menu  display 

CLEAR  @  DISP  USING  8600  @  DISP  USING  8650 

INPUT  WO 

IF  W0>9  THEN  Z0=W0  S  W0=0  ELSE  I4=W0 

IF  14=0  THEN  2630  !  return  to  main  menu 

ON  14  GOTO  2140,2200,2260,2310,2360,2400,2530,2590,2440 

1  5204  LI  A  settings (i 4=1) 

DISP  "  L I A  fullscale(mV)  jtiuie  constant  (ms)  ,  expand11; 

INPUT  V5,T5,G5@  V5=V5/  ( 1000*Q5) 

DISP  "delay:  tt  time  consts, extra  t  i  m  e  ( m  s 

PRINT  USING  8610  ;  1000*V5,B5, T5, 15, TS 

SOTO  2630 

!  water  depth, temp (i4=2) 

DISP  "column  depth  (mm) , temperature (C) "j 8 

D3=D1*D2/10A9 

PRINT  USING  3620  ;  Dl , D3 

PRINT  USING  S640  j  Wl 

SOTO  2630 

1  drive  voItages(i4=3) 

DISP  "drive  *or  P/T" ; V4; " / " ; V3; @  INPUT  M* 

IF  LEN(W$)>0  THEN  V4=VAL(WI)  @  V3"VAL (W$[P0S (W*, " 

PRINT  "P  drive:  "; V4; "7"  8  PRINT  "T  dr i ve; " ; V3 ; " V  " 

GOTO  2630 

1  comments  (14=4) 

DISP  "comments" 

INPUT  W* 

IF  LEN(W*)>Q  THEN  PRINT  "#";W*  @  GOTO  233J  ELSE  235C 

GOTO  2630 

1  various  constants(I4=5) 

CLEAR  S  DISP  "enter  new  values, then  CONTINUE"  Ci  PAUSi 

COPY 

G0r0  2630 


1  i 
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2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2510 


2540 
2550 
2560 
2570 
2580 
2590 
26  00 
2610 
2620 
2630 
26  40 


!  titl 
A*-A*[ 
DISP  " 
DISP  " 
1  tr  an 
A*=A*[ 
H$  =  "T 
DISP  " 
H$="H 
DISP  " 
PRINT 
Y3=Y3/ 
GOTO  2 
!  plot 
53=N0T 
IF  S3 
PLOTTE 
LIMIT 
GOTO  2 
1  posi 
DISP  " 
IF  LEN 
GOTO  2 
14  =  0  3 
1  end : 


e  (14=6) 

1,18] 

date";@  INPUT  WIS  A$  =  A$S<W$  @  GOSUB  8070 

run  #"  @  INPUT  W$@  A$=A$&W*  @  GOSUB  8070 

sducer  chars  o n 1 y ( i  4  =  9 ) 

1,64] 

is:"  S  DISP  H$;@  INPUT  W$@  A*=A$&H*&W*  @  GOSUB 

T  active  1 eng th (mm) , capaci tance (pF) " ; @  INPUT  D5 

is:"  S  DISP  H$;3  INPUT  W*@  A$  =  A*?.<HnW* 

H  active  length (mm)  ";@  INPUT  D6 

USING  8630  ;  A$,D5,Y3,D6 

10A12  @  D5=D5/1000  @  D6=D6/1000  !  SI  conversion 

630 

enabl e/di sabl e ( i  4  =  7) 

S3  @  13=0 
THEN  2560  ELSE  2630 
R  13  705 
3 

630 

tion (: 4=8) 

position-"; Z0;@  INPUT  W$ 
(W$) >0  THEN  Z0«VAUW#) 
630 

34=1  §  RETURN 

operator  input  procedure—- 


8070 
,Y3 


!2, 232, 30,  130  @  LOCATE  22,122,27,97  g  PEN  0 


59 


2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2300 
2810 

28  20 
2830 
2340 
2850 
2860 
2870 
2880 
2390 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 

29  3  0 
2990 
3  0  00 


!  begin:  peak  frequency  measurement  procedure  =  =  =  =  s:  =  =  = -•■=-  =  = 

PRINT  USING  8460  ;  A*  !  title 

CLEAR  @  DISP  "peak  frequency  measurement" 

DI5P  "estimate, search  r ange , stepsi ze" ; 2  INPUT  F7,WO,F9 

F4=IP(F7-W0) 

F5»CEIL(F7+W0) 

1  documentation 

PRINT  USING  8470  ;  TIME  ,  F7 , WO , F9 

1  warm up  reading 

IF  S5  THEN  V6  =  V3  ELSE  V6  =  V4 

F6^F7  @  T7  =  T8  @  GQSU8  6990  !  xmtr  on  ,long  delay 

GQSUB  7310  !  warmup  read 

IF  S3  THEN  59=1  @  GOSUB  5770  !  begin  plot  of  a  vs  f 

1  begin  : initialize  sampling  array 

F6=F4-3*F9  @  GDSUB  6990 

FOR  12=1  TO  7 

GOSUB  7310  !  read  amplitude  A2 

IF  12(7  THEN  F6=F6+F9  @  GOSUB  6990  !  xmtr  on  next  fraq  stei 

A(I2)=A2 

NEXT  12 

!  end  :  f  i  1 1  sample  arr ay 

1  begin: resonance  search 

S7,F1,F2,F7=0  !  peak  detect  flag  and  detection  storage 

FOR  F0*F4  TO  F5  STEP  F9 

1  fO  is  the  freq  at  the  center  of  the  sample  data  set 

F6=F0+4*F9  @  GOSUB  6990  !  xmtr  on  next  read  fraq 

1  curvature  at  f 0 

U2=(5*(A(7)+A(1))-3*(A(5)+A(3))-4*A(4))/B4 

DISP  USING  8480  ;  F0,U2 

IF  U2<0  THEN  2950  ELSE  3120 

1  slope  at  fO 

Ul-(S*(A(7)-A(l))+2*(A(6)-A(2>)MA(5)-A(3)))/29 

!   peak  frequency 

F3*F0-F9*U1/(2*U2) 

DISP  USING  3490  ;  F3 

IF    ABS(F3-F0XF9*R4    THEN    3010    ELSE    3120 


60 


3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3130 
3160 
3170 
3130 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
32E0 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
33  7  0 


E£ 


U0  = 

A7  = 

37  = 

IF 

IF 

PRI 

U9 

GOT 

!  2 

1  b 

S7 

1  e 

!  e 

DIS 

IF 

!  P 
FDR 
A(I 
NEX 
GOS 
NEX 
I  | 

IF 
IF 
IF 

1  < 
F7-= 
GOT 

!  t 

F7  = 

I  I 

PRI 
£4- 

RET 


ak  detected  at  center  of  sample  s.rr 

e a k  amplitude 

7*A(4)+6*(A(3)+A(5))+3*(A(2)+A(6))- 

0-UUUl/(4*U2) 

7+1  9  V7  =  V1 

5  THEN  V3  =  V6  ELSE  V4  =  V6 

7=1  THEN  F1=F3  ELSE  F2=F3 

T  USING  3500  ;  F3,1000*A7*Vi , 10Q0*U 

9  2  X0=FO 

3160 
d:peak  detected  block 
g  i  n  :  p  e  a  k  not  detected  block 

!  unset  peak  detect  -flag 
d:  peak  not  detected 
d:tests  at  FO 


ay 
2* 


A(7)+A 


\  i 


>/21 


2,1 


000*V1 


7<2  THEN  3180  ELSE  3240  !  end  \ 
ah  down  sampl ing  array 
10=1  TO  6 
)=A(I0+1) 

10 
B  7310  S  A<7)=A2 

FO 
d  search 

3  THEN  G03U8  5330  !  show  plot 
7<2  THEN  COPY  G  GOTD  3270  ELSE 
1-0  THEN  PRINT  "search  failed" 
e  good  reading 
1  3  F3^F9 

3330 
o  good  readings 

Fl  +  F2)/2  @  F8  :ABS(F2-F1)  !  avei 
n  a 1  results 
T  USING  3510  ;  F7,F3 

1  fast  refresh 
RN 
decanter  f r e q  m  m  t  procedurs  =  =  =  : 


e a  r  c  h  o  n  t  w  o  detects 


3310 
9  GO" 


a  q  e  o  ■ 


two 


9<?.d  1  112  3 
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3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 
3530 
3590 
3600 
3610 
3620 
3630 
3640 
3650 
3660 
3670 
3680 
3690 
3700 
3710 
3720 
3730 
3740 
3750 
7760 
3770 
3730 


!  beg 

CLEAR 

DISP 

Q7  =  0 

PRINT 

GQSUB 

A7*V8 

1  set 

F9  =  CE 

F4-MA 

F5-=F7 

DISP 

A4  =  A7 

DISP 

!'  bag 

F6-IP 

A8  =  A7 

1  beg 

F6=F6 

IF  F6 

GQSUB 

DISP 

1  che 

IF  A2 

1  uph 

PRINT 

BOTD 

1  d  o  w 

A8=A2 

IF  A2 

1  end 

beg 
IF  F6 
A  3  -  A  2 
GQSUB 
DISP 
IF  A2 
!  end 
F1=F6 
PRINT 

e  n  d 


in:  Q  measurement  procedure=~: 
3  DISP  "Q  measurement" 
"center  freq  =";F7 


ING  8570  ;  TIME,F7,Z0,3$ 
50  !  measure  peak  voltage 


US 

4950  !  measure  peak  voltage   v8   at  f7 
/Vl 

search    limits 
IKF7/R7) 
X(0,F7-20*F9) 
+20*F9 

"peak  : "; 1 0  0  0  *  A  7  *  V  1 
♦.707106731135  I  1/2  power  amplitude 
"half  power: "; 1000*A4*V1 
i  n : s  e  a  r  c  h  for  lower  1/2  point 
(F7) 


in  sloop  to  find  f6<fi 

-F9 

>F4    THEN    GQSUB    69*0    ELSE    4090    !    xfflti 

7310  !  measure  A2 
"q  ";  (A2-A4)#1000»V1 
ck  for  downhill  motion 
>A8  THEN  3620  ELSE  3650 
ill;  low  1/2  pt  obscured 

USING  8520  ;  F6,A2,A3  Q  COPY 
4110 
nhill- reset  check  v a 1 


0:1 


<  =  A 

lo 
i  n 

<F7 

i 

73 


lo 
-(A 

US 

SB 


4  THEN  3680  ELSE  3550 
op 
lo 
T 


i  Op! 

'HEM 


ve 
i 

;a2 

THE 
.br 

A4) 
G  3 
ch 


br 

F6 

las 

new 
•A4 
N  3 
ack 

/  (A 
530 
for 


t  1/2  point 

1  G  GQ3UB  6990  ELSE  4070 

mol  e 


i 


acke 
=  F6  + 
t  sa 
sam| 
)#1000#V1 
750  ELSE  36 
eted 
2-A8 

;  f 

Fl 


9n 


interpolate 
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3790 
3800 
3810 
3820 
3830 
3840 
3850 
3360 
3870 
3380 
3890 
3900 
3910 
3920 
3930 
3940 
3950 
3960 
3970 
3930 
3990 
4000 
4010 
4  0  2  0 
4030 
404  0 
4050 
4060 
4070 
4030 
4090 
4100 
4110 
4120 
4130 
4140 


1  b 
F6  = 

A3* 
!  b 
F6  = 
IF 
GOS 
DIS 
1  c 
IF 
!  u 
PR  I 
GOT 
1  d 
A8  = 
IF 

1  a 

!  b 

IF 

A8« 

GDS 

DIS 

IF 

1  e 

F2  = 

PR  I 

C7  = 

S4  = 

PRI 

34  = 

RET 

1  b 

PRI 

COP 

GOT 

1  e 


e  g  i  ,i : 
IP (F7) 
A7  !  c 
egi  n : 
F6  +  F9 
F6<F5 
UB  731 
P  "q  " 
heck  d 
A2>A3 
phill- 
NT  USI 
0  4110 

0  w  n  h  i  1 
A2 

A  2  <  =  A  4 
nd  loo 
egi  n  1 
F6>F7 
A2  !  s 
UE  731 
P  "  q  !  " 
A2>=A4 
nd  loo 
F6+(A2 
NT  USI 
F7/  (F2 

1  !  fa 
NT  USI 
1  !  fa 
URN 

a  g  i  n :  a 

NT  USI 

Y 

0  4080 

nd:  3  in 


search  for  hi  1/2  power  pt  F2 

heck  value 
loop  for  f 6>f 2 

THEN  GOSUB  6990  ELSE  4090 

0  !  measure  A2 

; (A2-A4)*1000*V1 

cwnhi 1 1  mot i  on 

THEN  3890  ELSE  3920 

hi  1/2  obscured 

NG  3520  ;  F6,A2,A8  @  COPY 

1-reset  check  val 


THEN  3950  ELSE  3320 


op. 

HEN 

ve 
i 

( A  2 
THE 
.br 

A4) 
G  3 
Fl) 
t  r 
G  8 


br 

F6 

1  as 

new 

-A4 

N  4 

ack 

/(A 

540 
i 

e  f  r 
550 
efr 


-  e 

6- 

sa 

am 

10 

0 

ed 

A3 

F 
al 
h 

0 
h 


t  1/2  point 

1  @  GOSUB  6990  ELSE  4090 

oipl  e 

pic?  aZ 

0  0  *  V  1 

ELSE  3960 

-interpolate 


i  t  y  factor 
7 


rror  processing 
NG  8560  ;  F6 


easurement  procedure32  —  --""""""""  — 
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4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4230 
4290 
4300 
4-3  1  0 
4320 
4330 
4340 
4350 
4360 
4370 
43S0 
4390 
4400 
4410 
4420 
4  430 
44  40 
4  450 
4  460 
4470 
4490 
4  490 
450  0 
4510 
4  520 
4  53  0 
761 
4540 
4550 


1  b 

CLE 

IF 

F6  = 

GOS 

IF 

DIS 

DIS 

PRI 

1  b 
FOR 
GOS 
Z4  = 

2  (I 
NEX 
1  e 
S7  = 

1  b 

2  4  = 
i 

U2  = 
Ul  = 
110  = 
U9  = 
12  = 
DIS 
IF 
1  h 
IF 
57  = 
Z3= 
1  c 

'  P 
FDR 
Z(I 
NEX 
GOS 
Z  (7 
IF 

i 

GO' 


egi 

AR 

S5 

F7 

UB 

33 

P  " 

?  " 

NT 

egi 

12 
UB 
ZO 
2!  = 
T  I 
nd: 
0 

egi 
Z(4 
par 
'(5* 
(3* 
(7* 
Z8 
Z4- 
P  U 
37 
it-' 
37  ' 
1  ! 
22 
on  t : 
ush 

0)  = 

T  I 

UB 

)=Z 

37 


n  :  w 
@  DI 
THEN 
@  T7 
7310 
THEN 
posi 
step 
USIN 
n :  i  n 
=  1  T 
7310 
@  GO 


avele 

SP  "m 

V6  =  V 

=  TB  @ 

!  wa 

S9  =  2 

t  i  on  = 

size" 

G  872 

i  t :  a  1 

0  7 


ngt 
in/ 
3  E 

GO 
rmu 

@ 

";Z 
;@ 
0  ; 

i  Z3 


h  measurement  prccedure=======: 

max  positioner" 


LSE  V6  =  V4 

SUB  6990  ! 

p  read 

GOSUB  5770 

0 

INPUT  Z9 
TIME,  20 ,  Z9 
sample  array 


select  T  or  P  drive 
xmtr  on, long  delay 

1  plot  start,  a  v s  z 


!  re 

SUB  7 

Z4  §  A ( I  2 


ad  amplitude  A2 

610  !  move  down  :9 
>=A2 


mm 


initialize  ar^ay 


bol 
A(7 
A(7 
( 4 ) 
XO 

a*u 

ING 

R  A 
doc 
HEN 

set 

sa 

nue 

dOW 

1  T 
(10 


earch  1 oop 

array  center 

i  c  fit 

)+A(i)  )-3*CA<5)+A(3)  )'-4*A(4)  )/84 

)-A(i)).+2*(A(6)-A(2))  +  (A(5)-A(3)  )  J/28 

+6*(A<S)+A(3) )+3*(A(6)+A(2) )-2*(A(7)+A 

=  24  S  A7»U0-Ul#U'i/(4*U2) 

l/'(2*U2)  !  predicted  optimum 


ins 


8660 


Z4,Z2,U1,U2 


BS(Z2-Z4XR5*Z9  THEN  4420  ELSE  4460 
u  (n  e  n  t 

end  search 

end 

estimate 
search 
n  sample  array 
0  6 
+1)  S  A(I0)=A(I0+1) 


4550  ! 

search 

v  e  peak 


310  !  read  a2(z0) 

0  A(7)=A2 
HEN  PRINT  USING  8730 


22,  1000*A7*Vi  ,Ui  ,U2    EL3C    G03UI 


0 


0  4320 

nd:  search  loop 
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4560 
4570 
4530 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4670 
4680 
4  6  90 
4700 
4710 
4720 
'30 


& 


4  7  40 

4  750 

4  7  6  0 

4770 

4730 

4  7  90 

4  3  00 

4810 

4320 

■1630 

4340 

4850 

4860 

19-2! 

4S70 

4880 

4890 

-1 9  )  0 

'"910 

4920 

4930 

49  4  0 


!  P 

PRI 

22  = 

PRI 

IF 

!  b 

IF 

IF 

DIS 

IF 

W0  = 

IF 

!  b 

19  = 

Hi, 

!  b 

CI- 

w  *- 

C3- 
C4  = 
C5- 
C6  = 

!  a 
!  b 
IF 
1  t 

Wi  = 
wo« 

C7  = 

T  p 

C3  = 

*  I  9 
1  d 
PRI 

COT 
!  e 
1  b 
PRI 

94r 


r i  nt  ^  i  na 1  est i  mate 

NT  USING  8730  f  Z2, 1000*A7*V1 ,U1 , 

(22+23)/2  !  average 

NT  "average*"; Z 2 ; "mm" 

S3  THEN  G03UB  5830  !  show  plot 

egin  data  revi  ew 

36  THEN  4630  ELSE  4640 

S5  THEN  V3=V6  ELSE  V4=V6  !  store 

P  CHR*  (11) j  "/4  count=";N9+2;@  INP 

LEN(W*)>0  THEN  4660  ELSE  N9=N9+2 

VAL(WI) 

WOO  THEN  N9  =  WQ  ELSE  4910 

eginzaccept  measurement 

19+1 

L9-22/1000  !  meters 

egin:  sums 

C1+N9  !  -:< 

C2+N9*N9  ! 

C  3  +  L  9  !  *  y 

C4+N9+L9  ! 

C5+L9*L9  ! 

I9*C2-C1*C1  !  del  factor 

n d :  sums 

egimcompute  sound  s p d 

C6/0  THEN  4810  ELSE  4370 

wc  points  or  more 

(I?*C4-C1*C3)/C6  !  slope 

(C2*C3-C1*C4)/C6  '  intercept 

4»W1*F7 

I9>2  THEN  4860  ELSE  4870 


L2 


ADC  r 
UT  W* 

§  GGT 


esui  t s 

G  46  80 


a  n 


y 


^y 


o + w  i* n 


■  2  *  W  0  *  i 


18,  400  )*WJ 


4 * 3 D 9.  (  (C5  +  I9#W0*I 
/C6)*F7 

o  c  u  m  e  n  t 

NT  USING  B690  ;  19. 

0  4930 

nd  :  accept  measuremt 
egin:  discard  measuremt 
NT  "  m  m  t  discarded" 

1  S  RETURN 

nd:  wavelength  measurement  p  r  o  c :?  i 


■  W 1  *  C  4  +  2  *  W  0  * .; 
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4950 
4960 
4970 
4990 
4990 
5000 
5010 
5020 
5030 
5040 
5050 
5060 
5070 
5080 
5090 
5100 
5110 
3120 
5130 
5140 
5150 
5160 
5170 
5 1  3  0 


!  b 
IF 
F6  = 
GOS 
1  b 

so= 

F6  = 
GOS 
F6  = 
VS  = 
GOS 
F6^ 
V8  = 
GOS 
F6  = 
V3  = 
GOS 
F6  = 
VB» 
GOS 
V8  = 
S6== 
RET 
1  s 


eg  l 

S5 

F7 

UB 

egi 

S6 

F7- 

UB 

F6  + 

-(3 

UB 

F6  + 

V8  + 

UB 

F6  + 

V8  + 

UB 

"6  + 

ve  + 

U3 
V8- 

so 

URN 
nd; 


n:  5 

THEN 

@  T7 

7310 

nifi 

@  S6 

2  @ 

7310 

1  @ 

*V27 

7310 

i  a 

12*V 

7310 

1  @ 

17*V 

7310 

1  S 

12*V 

7310 

3*V2 


peak  voltage   measurement  prcceoure: 
=  V3  ELSE  V6  =  V4  !  select  T  or  P  drive 

@  GOSUB  6990  !  :<  m  t  r  on  ,  1  o n g  delay 
warmup  read 

point  mint  of  peak  amplitude  v2 
!  ADC  olf 

UB  6990  !  set  xmtr  Fn-2 
read  v2 (Fn-2) ;  1st  term 
UB  6990  !  xattr  on  Fn-1 

!  first  term 
read  at  Fn-1 
UB  6990  !  xmtr  on  Fn 
5  !  2d  term 
read  at  F n 

UB  6990  !  xmtr  on  Fn+i 
5  '  center  term 
i  read  at  Fn+1 
GOSUB  6990  !  xmtr  on  Fn+2 
2/35  !  4th  term 
!  read  at  Fn  <■  2 
/35  !  result  at  f7 
restore  ADC 


-pt 

V6 

=  T8 
i 

ve- 

=  0 

GOS 
i 

GOS 

35) 

i 

GOS 

2/3 
i 

GOS 
2/3 


peak  voltage 


mint1 
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5190 
5200 
5210 
5220 
5230 
5240 
5250 
5260 
5270 
5230 
5290 
5300 
5310 
5320 
5330 
5340 
5350 
5360 
5370 
5380 
5390 
5400 
5410 
5420 
5  430 
5440 
5450 
5460 
5470 
34  90 
3490 
5500 
3510 
3320 
3530 
3540 

si  U  J  U 

J  J  o  u 
5S7A 


!  begin:  calibrate  hydrophone  procedurs  =  =  : 
W1=C7/F7  !  wavelength 


H1=FNH(PI*D5/W1)  !  len 
H2=FNH(PI*D6/W1)  !  len 
PRINT  USING  8700  ;  A*, 
!  begin; preamp  cal i  bra 
CLEAR  @  DISP  "calibrat 
13=4  §  GOSUB  8120  !  sm 
DI3P  "preamp  calibrate 
F6  =  F7  @  V6  =  V5*S8/G2.  @ 
GOSUB  7310  !  read  LIA 
G2*G7*A2 

PRINT  USING  3740  ;  FNL 
!  end :  preamp  cal i  brat 
1  T'H  measurements 
18=3  B  GOSUB  8120  !  sw 
DISP  "TJH"  @  INPUT  W$ 
GOSUB  4950  !  5-pt  mmt 
E3=V8*H2/G2  !  Eth 
J3=V1*2*PI*F7*Y3  !  amp 
PRINT  USING  3750  ;  VI, 
1  end  T:H 
1  b  9  g  i  n  P  }  H 
I8=-2  3  GOSUB  3120 

Disp  "p:h"  e  input  u$ 

1  P  dr  i  vs  adjust 
F6^F7  @  V6--V4  3  T7  =  T3 
GOSUB  7310  !  warmup  re 
'  begin:  adjust  loop 
GOSUB  6990  !  xmtr  on 
GQSUD  7310  @  E2=V2*H2/ 
DISP  USING  37 10  ;  1000 
IF  ABS(E2-E3) >R3»E3  TH 
1  end   adjust  loop 
V4  =  V6 
PRINT 
!  end 


gth  aysrage  corrn 

gth  average  corrn 

FNL(Hl)  ,FNL(H2) 

ti  on 

ion" 

itch 

11  @  WAIT  T 3 * 3 

T7=T3  S  GOSUB  699( 

cut 

(G2) ,1000*V1 
i  on 

itch 

of  v8 

eras 

1  0  0  0  *  J  3  ,  1 0  0  0  0  0  0  *  e  ; 


for 
for 


mtr 


Oil 


(3  GOSUB  6990 
a  d  i  n  g 


;<  m t  r  on 


G  2  '  read    v 2 ,  correct 
*V6, 1000*V2/62, 1000000; 
EN  V6*FNR(V6*E3/E2)  @  I 


p  r 


Q'r3 


&  a  m  p 


"P  drive: ";V4; "V 
:  P  drive  adjust 

GOSUB  4950  !  vS  mmt 

E2=V8*H2/G2  !  Eph 

PRINT  "  Eph="; 1000000*1 

!  end  P)H 
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5600  !  begin  P>T 

5610  DISP  "PH"  3  I3  =  -l  3  GOSUB  8120 

5620  GOSUB    4950    !    v8    mint 

5630  E1=V8*H1/G2  !  Ept 

5640  PRINT  "  Ept  =  " ;  1 000000*E1 ; "uV" 

5650  !    end    P}T    mints 

5660  !  begin  sensitivity  calc 

5670  J7=PI»D3*F7/(K0*C7*C7*Q7) 

5630  PRINT  USING  8760  ;  F7,C7,Q7 

5690  PRINT  USING  8770  ;  D3,K0,J7 

5700  Mi-8QR(J7*E1*E3/(E2*J3)) 

5710  M2=SQR(J7*E2*E3/(E1*J3) ) 

5720  PRINT  USING  87B0  ;  FNM  (  Ml )  ,  FNM  <  M! 

5730  '  Lsnd  sensitivity  comp 

5740  54=1  @  RETURN 

5750  '  end:calibrate  phonesssas======: 


switch 
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5760 
5770 

5780 
5790 
5800 
5810 
5820 
5830 
5840 
5850 
5860 
5870 
5830 
5690 
5900 
5910 
5920 
5930 
5940 
5950 
5960 
5970 
5980 
59  90 
6000 
6010 
6020 
6030 
6040 
6050 
6060 
6070 
6080 
6090 
6100 
6110 
6120 
6  1  3  0 
6140 
6150 
6160 
6170 
6180 
6190 
6200 
6210 


:  =  =  =  =  PL0T  PRQCS  =  -=«  =  = 

n:  plot  initialize  procedure5 
load  plotter";®  INPUT  W$  !  du 
l(W$) >0  THEN  33=0  @  RETURN  !  a 
@  FRAME  @  PENUP 


n 
n 

9 

0 
>B 
<B1 


pi  o 
:  s 

:  s 

999 
1  T 


2  T 

T 
5  T 
4  T 


t  in 
how 
cali 
9  2 
0  13 
Wl  = 
HEN 
HEN 
HEN 
HEN 


i tial ize  procedure; 
plot  procedure"  — 

ng 

B2,B5=-B1 

Y(I0) 
B2«W0 
B1=W0 
B3-W1 

B4  =  Wt 


!  beg 
DISP 
IF  LE 
PEN  1 
RETUR 

end 

beg 

beg 
Bi,B4 
FOR  I 
W0  =  X( 
IF  WO 
IF  WO 
IF  HI 
IF  Wl 
NEXT 
DISP 
DISP 
!  X  a 
B3  =  5 
WO=5* 
B8=(C 
IF  B8 
IF  B8 
1  too 
K  3  =  1 

B3=B3*B(K3)  @  K3=FNK(K3) 
W0=5*83 
B  8  ••=  (  C 
IF  B3 
!  end 
1  beg 
K3  =  3 
B  3  ■-  B  3 
IF  B3 
W0  =  5* 
BB»(C 
IF  BS 
!  and 
1  ?nd 
3 1  =  I P 
B2  =  B1 
GOTO 
1  met 
B3=(B 
BB  =  10 


(nmy 
bort 


i  nput 
plot 


a  t  a  range  x :  " ;  8 1 ;  B  2 
y:  ";B4;B5 
s  scale 
x-interval 

!  major  divisions 
L  ( B  2 / W 0 ) -  I P ( B 1 / W 0 ) )*5 
0  THEN  6010  ELSE  6000 
0  THEN  6080  ELSE  6160 
any  divisions 


EIL(B2/W0)-IP(B1/W0) )#S 

>50  THEN  6030  ELSE  6160 
too  many  divisions 
in  too  few  divisions 

/B(K3)  @  K3«FNK(K3) 

>0  THEN  6120  ELSE  6200 
B3 

EIL(B2/W0)-IP(B1/W0) )*5 
<20  THEN  6100  ELSE  6160 

too  -few  divisions 

if 

CB1/W0)  *W0  !  low  :<  bdry 
+  B 3 * B 8  !  hi  ;<  bdry 
6250 

hod  failed-set  de-fault 
2-B1)  /10 


:<  scaling 
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6230  DISP  "de 

!  end  de 

I  end  x 

!  begin 

B6=. 002 

W2=5*B6 


6240 
6250 
6260 
6270 
6280 
6290 
6300 
6310 
6320 
6330 
6340 
6350 
6360 
6370 
6380 
6390 
6  400 
6410 
6420 
6  430 
6440 
6450 
6460 
6470 
6480 
6490 
6500 
6510 
6520 
6530 
6540 
6550 
6560 
6570 
6580 
6590 
6  6  0  0 
6610 
6620 
6630 
66  40 
6650 
6660 
6670 
6680 
6690 
6700 


B9=  (CEIL 
IF  B9>40 
IF  B 9 <  1 0 
1  begin 
K3  =  3 

B6=B6*B( 
B9= (CEIL 
IF  B9>40 
'  end  to 
!  begin 
K3=l 

B6=B6/3( 
IF  86 >0 
B?=  (CEIL 
IF  B9<  10 
1  end  to 
1  end  if 
B4-1P (B4 
B5=B4+B6 
GOTO  654 
1  method 
86- (B3-B 
33=10 
DISP  "de 
end  d  e 
!  end  y  - 
SCALE  Bl 
1  end  sc 
!  plot  p 
OUTPUT  7 
K3=I3-7 
IF  K3<1 
FOR  10=1 
PLOT  X(I 
NEXT  10 
OUTPUT  7 
K3=K3+1 
FOR  I0=K 
PLOT  X(I 
NEXT  10 
OUTPUT  7 
CSIZE  4 


■fault  x -seal ing" 
fault  scaling 
seal ing 
y  scaling 
1  mV 

(B5/W2) -IP  CB4/W2) )*5 

THEN  6320  ELSE  6310 
THEN  6380  ELSE  6450 
too  many  divisions 

K3)  @  W2=5**B6  a  K3«FNK(K3) 
(B5/W2) -IP  (B4/W2) )*5 

THEN  6340  ELSE  6450 
o  many  divisions 
too  few   divisions 

K3)  9  W2=5*BA  @  K3=FNK(K3) 
THEN  6420  ELSE  6490 
(B5/W2) -IP (B4/W2) ) *5 

THEN  6400  ELSE  6450 
o  few   divisions 

/W2) *W2 

*39  !  hi  y  bdry 

0 

' f a  i 1 -  s  e  t  default  y - s c  a  1 i  n  g 

4  )  / 1  0 

■fault  /-scaling" 

fault  y  scaling 
seal ing 
,B2,B4,85 
al  i  ng 
c  i  n  t  s 

05  ;  "SMo;  " 

THEN  K3=0  @  GOTO  6640 

TO  K3 
0  >  ,  Y  ( 1 0 )  ,  0 


;k* 


3  TO  13 
0  )  ,  Y  (  I  0  )  ,  0 

05  ; "SM;  " 
S  FXD  1,3 
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6710  LAXES  B3,B6,B1,B4,5,5 

6720  N4  =  N4+ 1  @  H$=A*C33 , 371  @  PRINT  "plot  " ; H * j " # " ; N 4 

6730  SETGU 

6740  MOVE  0,62  @  LDIR  90  @  LQRG  6 

6750  LABEL  USING  "K"  ;  " amp  1 i tude (mV  rms)" 

6760  LDIR  0  @  LORG  4  @  MOVE  72,12 

6770  IF  S9=2  THEN  W*  =  " dep th ( mm) "  ELSE  W*="f requency (Hz ) " 

6780  LABEL  USING  "K"  ;  W$ 

6790  MOVE  72,5  (3  CSIZE  5 

6B00  IF  59-2  THEN  LABEL  USING  3670  ;  F6,SI,H*,N4  ELSE  LABEL  USING 

8630  ;  Z0,S*,H*,N4 
6310  SETUU  9  PENUP 
6820  !  end:  point  plot 
6830  !  begin:  curve  fit  plot 
6840  B3=(B2~B1)/100  ■ 
6950  FOR  W1=B1  TO  B2  STEP  B3 
6860  W2=W1-X0 

6870  W3-i000#Vi#(U0+Ui#W2/U9+U2*(W2#W2) /(U9*U9) ) 
6880  IF  W3>B4  AND  W3(B5  THEN  PLOT  W 1 , W 3 
6390  NEXT  Wl 
6900  PENUP 

6910  MOVE  0,0  9  PEN  0 
6920  !  end:  curve  tit  plot 
6930  !  begin:  clear  plot 
6940  13=0  6  S9=0 
6950  \    end:    clear   plot 
6960  RETURN 
6970  !  end:   plot  procedure  =  =  =  =  =  =  -  =  =  —  -  =  --  =  =  =  =  -  -  =  =  -  =  s  =s  =  =  =  = 
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6930 
6990 
7000 
able 
7010 
7020 
7030 
7040 
7050 
7060 
7070 
7030 
7090 
7100 
7110 
7120 
7130 
7140 
7150 
7160 
7170 
7130 
7190 
7200 
7210 
7220 
7230 
7240 
7250 
7260 
7270 
7200 
7290 
7300 


begin 
input 
si  oft 
IF  S6  T 
1  begin 
!  autom 
OUTPUT 
31  =  0  @ 
1  begin 
WAIT  2* 
ENTER  7 
IF  A5>. 
1  LI  A.  o 
IF  ASM 
GOTO  71 
IF  A5<. 
!  LIA  o 
V  6  =  M  I  N  ( 
GOTO  71 
1  reset 
OUTPUT 
OFF  TIM 
ON  TIME 
PRINT  3 
1  hold 
IF  31  T 
1  end  : 
'  begin 
OUTPUT 
SI  =0  <3 
1  e  n  d  :  n 
T7-0  @ 


=  MEASUREMENT  PROCEDURES*""""™-"""™™"-- 
:  amplitude  read  setup  pr Qcedure-"=""=""""  =  = 

:  3  6  ,  F  6  ,  V  6      out:    tinier    #1    set    for    read    enable,    read    sn 


HEN  702 

:  ADC  b 

at i  c  dr 

707  USI 

ON  TIME 

:  1  ocp 

T5 

09  ;  A5 

8  THEN 

ut  too 

. 12  THE 

70 

1  THEN 

ut  too 

10,FNR( 

70 

HP3314 
707  USI 
ERtt  1  3 
Rtt  1,T8 
$[1,1]; 
for  r  e  a 


0  ELSE  7250 

lock 

i  v  e  adjust 

NG  3800  ;  V6,F6  !  HP3314  drive 

Rtt  1,T7+I5*T5  G0SU8  7550  I  rea 

to  adjust  V6  until  LIA  output 


on 
d  di  sabi  s 
is  in  limits 


1  read  LIA  magnitude-ignore  read  d i s a b 1 
7100  ELSE  7130 
hi  gh 
N  V6=FNR(V6*.5>  ELSE  V6-FNR(V6< 


7150  ELSE  7290 
1  ow 
76*15) ) 


and  retast 
NG  8810  ;  76 

S 1  =  0 
+(I5-2)*T5  GOSUB  7550 
"  drive"; V6; "V" 
d  enable 
0  ELSE  .7220 
ck 

C  block 
NG  3300  ;  V6,F6 


i*. 7) 


here's  the  loop  exit  point 


ADC  b  1  0! 
:  no  AD1 

707  USI! 

ON  TIMERS  i,T7+I5*T3  GOSUB  7550  !  read  disabl 

o  ADC  block 

RETURN 


end:  amplitude  read  setup  procedure  ::  : 
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7310 

7320 

7330 

7340 

7350 

7360 

7370 

7390. 

7390 

7  400 

7410 

7420 

7430 

7440 

7450 

74  60 
7  470 
7490 
7490 
7500 
7510 

75  20 
7530 
7540 
7530 
7560 
7370 
7590 
7390 
7  6  0  0 


begi 
in:  9 
hold 
IF  91 
!  read 
TRIGGE 
ENTER 
V2=A5* 
ENTER 
A2=V2/ 
!  chec 
IF  A5> 
1  al  a.r 
PRINT 
!  snd : 
U 


n  :  a  m  p 
1  ,V0 

for 
THEN 

enab 
R  723 
709  ; 
V5  ! 
723  ; 
VI 

k  for 
1  THE 
m  -ov 
"LI  A 

over 
SING 
AND  S 


litude  read  procedure-  =  =  =  "  =  =  -  =  =  : 

out:  v2,a2 

read  enable 

7350  ELSE  7330 

le 

!  measure  ;<  m  t  r  volt 

A5  !  LIA  magnitude 
volts  recei  ver   out 

VI  !  xmtr  voltage 


aga 
outp 

put 


ut 


DI 

IF  S3 

!  graphics 

13=13+1 

IF  S9--2  THE 

Y( 13) =1000* 

1  end  graph 

RETURN 

end  :  amp 
begin: r  e a 
called  by 

OFF  TIMER* 

31-1 

RETURN 

1  end:  read 


LIA  overload 
N  7430  ELSE  74 
er 1 oad 

OVERLOAD"  9  PR 
load  test 
9590  ;  F6,Z0,i 
9  AND  I3<50  TH 
storage 


50 

INT  F6;A5;V1  S  BEEP  55,2000 


000*V1 

EN  743 


,100 
0  EL 


o*v 

SE 


2 
7530 


N  X  ( 13) -Z0  ELSE  X (I3)=F6 

V2  !  mV 

i  cs  storage 


litude  read  procedun 

d  enable  procedure-- 

timerSl-resets  read 


enable  f 1  a g  s 1 


enable  procedure==~=- 
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7610 
7620 
7630 

7640 
7650 
0"  @ 
7660 
7670 
7690 
7690 
7700 
7710 
7720 
7730 
774() 
7750 
7760 
7770 
7780 
7790 
7S00 
7810 
7820 
7830 
7340 
7350 
7360 
7870 
7330 
7890 
79  00 
7910 
7920 
7930 
7940 
7950 
7960 
7970 
7930 
7990 
2  0)0 
3010 
3020 
3030 
8  0  40 
8050 
8060 


!  bsgin:  move  manipulator  down  Z9  mm"  =  =  =  =  =  =  =  =  =  =  "-  =  =  =  =  -3~  = 

!  input:  :9 

ENTER  720  ;  NO  !  counter  ref 

IF  ABS(Z9K.5  THEN  7630 

IF  Z9>0  THEN  OUTPUT  709  ;"CLS1"  @  N3=l  ELSE  OUTPUT  709  ;"CLS 

N3=-l  3  Z9=-(Z9/R8) 

ON  TIMER#  2 ,  IP  ( Z9/R2-T2)  GOSUB  7700  !  set  stop  time 

S1,S2=0  !  read  off,inmotion 

RETURN 

end:  move  manipulator  procedure"""3 -:"  =  =  """"-::::;::::::=:-~ 

bsgin:  manipulator  stop  procedure3"  =  """""~":-""u;" 

input  zQ,z9,nO,n3 

output  revised  z 0 , s t e p  z  8 
N 3 > 0  THEN  OUTPUT  709  ;"QPN1"  ELSE  OUTPUT  709  ;"QPNO" 

F  TIMERS  2 

TER  720  ;  Ml 

1 oop  unt i 1  ctr  stops 

TER  720  ;  N2 
N2-N1  THEN  7790  ELSE  N1=N2  3  GOTO  7760 

motion  stop 

=N3*(N1-N0) /RO 

=Z0+Z8 

3P  USING  3790  ;  Z0,Z9 
S2-0  !  flag  notion  stop 
!  auto  drive  control 
IF  36  THEN  7860  ELSE  7950 
!  beginssdjust  drive 
ENTER  709  ;  A5 

IF  A5>.95  THEN  7890  ELSE  7950 
!  LIA  out  too  high 

ASM.  12  THEN  v"6  =  FNR  ( V6*.  5)  ELSE  V6  =  FNR  (V6*. 75 


IF 
OF 
EN 


IF 

i 

Z3 

zo 

DI 


I 

OUTPUT  707  USING  8810 

PRINT  "drive: ";V6; "V" 

WAIT  I5*T5*T3 

GO 

ON 

T7 

RE 


V6 


TO  7  3  70 
TIMER*!  i,T7  +  I5*T5  GOSUB  7530  i  set  delay  to  r 

=  0 

TURN 
1  e n d  :  m a n  i  p u  1  a t o r  stop  procedure"3-3""""3"3 

begin:  manipulator  relocate  procedure3"33"33 
CLEAR 

IF  10=0  THEN  DI3P  "position";^  INPUT  ZO 
DI3P  "move  to" 5@  INPUT  Z4 
Z9  =  Z4-Z0  !  in  m 
GOSUB  7610  !  move  z9  mm 
34=1  @  RETURN 
'  end:  manipulator  relocate  procedure333-33333333 


sad  enabl 
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8070 
8080 
8090 
8100 

8110 
8120 
8130 
8140 

8150 
8160 
8170 
8130 
8190 
8200 
8210 
8220 
8230 
3240 
8230 
3260 
3270 
0230 
8290 
3300 
9310 
9  TIC 
8330 
3340 
3350 
3360 
8370 
3330 
3  3  9  0 
8400 
3410 
3420 
2430 
3440 


!  begin:  title  string  jListifier  =  =  =  =  -  =  =  =  =  =  -  =  =  =  =  =  =  =:-  =  =  ---- 

I0=LEN(A$) 

IF  10  MOD  3 2 > 0  THEN  A**A*&"  "  9  GOTO  3030 

RETURN 

!  end  title  ju,5ti-fier~  =  =  =  -  =  -  =  :::  =  =  =  "  =  =  =  :=  =  :=  =  ~~:==:-~  =  =  "::i:=--=  =  "; 

!  ADC/xmtr  toggle  routine"  =  =  ~  =  =  =  =  =  =  =  =  =  =  "  =  -  =  r" ===  =  =  =  =  =  ■ 

IF  13=0  THEN  3140  ELSE  3170 

1  operator  call 

S4  =  l  !  -fast  re-fresh  menu 

CLEAR  9  DISP  "1  =  P}T.  2-PJH  3-T}H  4=cal";@  INPUT  13 

!  process  command 

IF  1 8  >0  THEN  S6=l  ELSE  S6  =  0  @  1 3  =  - 1 3  !  I3>0  turns  ADC  on 

18-13  MOD  5 

IF  18=0  THEN  3420 

OUTPUT  707  ; "APOVO" 

ON  IS  SOTO  3230,3290,3330,3330 

1  P>T 

S5=0  9  S*="P>T" 

OUTPUT  706  USING  8820  ;  "B1A1" 

OUTPUT  714  USING  3320  ;  "32A2" 

GOTO  8420 

1  P)H 

55=0  @  S*="P}H" 

OUTPUT  706  USING  3820  ;  " B 2  A  I " 

OUTPUT  714  USING  3820  ;  "B2A2" 

GOTO  3420 

1  t:h 

35=1  @  S$*"T>H" 

OUTPUT  706  USING  3820  ;  "B2A2"  ! 

OUTPUT  714  USING  3320  ;  "B3A3" 

GOTO  8420 

!  p  r  e  a  :ji  p  calibration 

S5=0  @  S*="P]cal " 

OUTPUT  706  USING  8320  ;  "BJA1" 

OUTPUT  714  USING  3820  ;  "B7A1" 

PRINT  USING  "  '}•'/"  ;  5* 

18=0  @  RETURN 

!  end:  ;<  mtr/adc  toggle  r  out  i  n £==========  =  =  =  =  =  =  -.izm  ■.-.■=  *■.:-.--  --.-. 
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8450  !  Images  for  formatted  outputssss=a==!S  =  as  =  =  ====  =  :I!S  =  =  =  =  =  sa-  = 

6460  IMAGE  /////K 

8470  IMAGE  //"  t  i  me:  ",  60/  "center  freq  measur  ement "  /  "est  i  mat  e  :  l!  ,  50/ 

"range: " , 30/ "stepsi za : " ,30/ 

3430  IMAGE  »,"*", 4D , 2X , S . 4DE 

8490  IMAGE  #,90.20 

3500  IMAGE  /  "peak  :  "  ,  50.  20  ,  "  Hz  "  /5X3D.  50  ,  "  mV"  /5XS  .  3DE  ,  "  mv7HzA2"  /  ">:m 

tr:  "  ,K, "mV" 

3510    IMAGE    /"resonance    center    f requency'75D. 2D, "+/-" ,20. 3D , "Hz" 

8520  IMAGE  /"1/2  power  pt  obscur ed " / 40 , 2 ( XK ) 

3530  IMAGE  "low  1/2:",5D.D 

8540  IMAGE  "high  1/2: "  ,50. D 

3550  IMAGE  "resonance  quality  factor :", 2D. 2D 

8560  IMAGE  /"search  limit  exceeded" ,50 

8570  IMAGE  / "  t i me:  "  ,  60/  "Q  measur ement "  / " center  f req j " ,5D. 2D/"dept 

h  "  ,  4  D .  0  /  K  / 

3580  IMAGE  "  o"  ,  4D  ,  X4D . D , 60. 20 , 3D. 50 

3590  IMAGE  /K/  "Fn  =  " ,  K, "  +  /-",  2D.  2D/"C'n  =  "  ,K, "+/-",  2D.  20/  "Qn  =  -•  "' f-: 

T )  =» " ,  K , "  d  B "  / "  M  ( H )  - " ,  K , "  d  B " 

3600    IMAGE    /"Orsturn    to    main    i=LIA"/"2=watsr  3    drive"/" 

4«comments  5=constarits" 

3610    IMAGE    "lockin    analyzer "/"input    " , X2D. 4D, "mV" /"expand" , X3D/"t 

line    ",X5D,"x",2D,"+",X5D,"ms" 

56  20  IMAGE  "depth  (mm) : " ,K/"volume (mA3) : ",0.60 

3630  IMAGE  K/"T  1 ength=" , K , "mm" / "capaci tance=" , K , "pF" / "H 

,K,  "mm"/ 

8640  IMAGE  "water  temper ature " /2D. 0 , "  C" 

2650  IMAGE  "6=titla         7=plot  on/off  "/ "8  =  posi ti on 

d  c  r  specs" 

3660  IMAGE  "  A "  ,  40 . 20 , 2 X K / "  "  ,2 (XS.3DE) 

E670  IMAGE  40, "H:  "  ,K,7X ,K, "#" ,20 

3680  IMAGE  30.0,  "mm  "  ,K,6X,K, "#" ,20 

B690  IMAGE  /2D,"  pts"/  "spd=" ,  40. 2D, "  +  /-" ,20. 20, "m/s"/ 

3700  IMAGE  K/"T  corrn: " , S2D. 20 , "dB" /"H  corrn: " ,S2D.2D,  dB" 

3710  IMAGE  "pda",3D.2D,2<X.4DE) 

3720  IMAGE  /  *' 1 1  me :  ' ,  60  , "  max/min  1  ocator  "/"start :",  4D.  D, "  step:", 

30.0/ 

2730  IMAGE  "mi  n  /max :  "  ,  40 .  D  ,  " mm "  / 3X30  .  50  ,  "  n  .; "  ' si  ops:  " , S.  4DE/ ,!cur\ 

:  "  ,  S  .  3  D  E  / 

374C  :MAGE  "preamp  gai  n : " ,  3D.  20,  "dB   @  ",K,"mV" 

5750  IMAGE  "  V t  =  " , K  , " V" / "  1 1  =  "  ,K , "mA"/ "  Eth - " , K  , " uV " 

6760  IMAGE  // "f r eq :  " , 50. 20/ "sound  spd:  " ,  40. 2D/"Q: " ,3D. 3D 

B770  IMAGE  "volume:  "  ,K/ "densi ty : " , 40. 20/ " Jn : " ,K/ 

8780  IMAGE  "T: " , 40. 2D , "dB  re  i V/uPa" / "H: " , 4D. 20, "dB  re  IV/uPa"// 

8790  IMAGE  "stop: " ,40.20,50.40 

8800  IMAGE  " AP" , 20. 50 , " VQ  PR",  40, "HZ" 

8810  IMAGE  "AF",2D.5D, "VQ" 


!  n  q  t .": ' 


,■   ,  "mm" 
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8320 
9330 
8840 
3850 
8860 
8870 


IMAGE  »,4A 

I  end:  i  mages=======: 

!  begin:  end  program 
OUTPUT  709  ; "RS" 
OUTPUT  707  ; "APOVO" 
END 


routine  —  — 
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